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^sj , We prove an asymptotic stability result for the water wave equations linearized around 

small solitary waves. The equations we consider govern irrotational flow of a fluid with constant 
density bounded below by a rigid horizontal bottom and above by a free surface under the 
influence of gravity neglecting surface tension. For sufficiently small amplitude waves, with 
waveform well-approximated by the well-known sech-squared shape of the KdV soliton, solutions 
• of the linearized equations decay at an exponential rate in an energy norm with exponential 

, weight translated with the wave profile. This holds for all solutions with no component in (i.e., 

symplectically orthogonal to) the two-dimensional neutral-mode space arising from infinitesimal 
translational and wave-speed variation of solitary waves. We also obtain spectral stability in an 
unweighted energy norm. 
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Asymptotic linear stability of solitary water waves 



1 Introduction 

The discovery of solitary water waves by J. Scott Russell in 1834 was a seminal event in nonlinear 
science. Russell's observations gave him immediate confidence in the significance of these waves, 
and led him to carry out an extensive program of experiments investigating solitary waves and their 
interactions |40j. But mathematical understanding was slow to develop. The first significant steps 
forward were made by Boussinesq [HI [9l [1^ and Rayleigh |38] by carefully balancing long-wave 
and small-amplitude approximations. The simplest useful model (derived by Boussinesq already 
in 1872, see [TOl p. 360] and [30j) is the famous Korteweg-de Vries equation [25j. Its sech^ soliton 
solution approximates the shape of small-amplitude solitary water waves. 

Given the status of the KdV equation as an approximate model, it is important to understand 
whether the soliton solutions of the KdV equation are approximations of some solutions of a more 
exact water wave model with similar properties. In this paper, we focus on questions of stability 
for exact solitary wave solutions of the Euler equations that govern incompressible and irrotational 
motions of an inviscid, constant-density fluid of flnite depth. The fluid occupies a two-dimensional 
domain whose lower boundary is a flat rigid bottom and whose upper boundary is a free surface 
that forms an interface with air of negligible density and viscosity. Surface tension on the free 
surface is neglected. 

For these water wave equations, the existence of solitary wave solutions with shape well- 
approximated by the KdV soliton was proved by Lavrent'ev [57], Friedrichs and Hyers [T3] and 
Beale [IJ. If the surface tension is positive and small, finite-energy, single- hump solitary waves 
are not known to exist, and indeed, exact traveling waves approximated by the KdV soliton may 
not exist without 'ripples at infinity' [21 |32]. For large surface tension, solitary water waves of 
depression exist j42| . but the relevant physical regime corresponds to water depth less than 0.5 cm. 

Explaining the stability of solitary water waves mathematically remains a very challenging prob- 
lem, despite considerable physical and numerical evidence. Remarkably, a valuable step forward 
was made already by Boussinesq [9l [10] , who argued for their stability based on a quantity he called 
the 'moment of instability,' which he showed was invariant in time based on the KdV approxima- 
tion. Over a century later, Benjamin [3] made use of the same quantity as a Hamiltonian energy, 
constrained by a time-invariant momentum functional, to develop a rigorous variational method 
to prove orbital stability for the set of solitary- wave solutions of the KdV equation. Benjamin's 
arguments were improved and perfected by Bona [5|. 

Variational methods for orbital stability and instability in Hamiltonian wave equations, based 
on the use of energy-momentum functionals, were subsequently greatly advanced by many authors. 
Notably, the general theory of Grillakis et al. [20[ [2Tj has been applied extensively to many physical 
systems. Using variational methods of this type for the case of solitary water waves of depression 
for the Euler equations with large surface tension, orbital stability conditional on global existence 
was obtained by Mielke [55] and Buffoni [llj. For small surface tension, such variational stability 
results have also been obtained recently by Groves and Wahlen |22] for oscillatory traveling wave 
packets of finite energy (also called solitary waves by several authors). 

For solitary waves with zero surface tension, however, it appears hopeless to study stability 
using variational methods based on constrained minimization. As remarked by Bona and Sachs [6], 
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the usual energy-momentum functional is highly indefinite in this case — The second variation lacks 
the finite-dimensional indefiniteness property key to the success of current variational methods. 
Regarding the stability of solitary waves with zero surface tension, the only existing rigorous work 
appears to be the recent paper of Lin ^28j, which addresses the linear instability of large waves close 
to the wave of maximum height. 

The present study involves a direct analysis of the Euler equations linearized about a small- 
amplitude solitary wave solution. The linearized equations have a natural two-dimensional space 
of neutral modes arising from infinitesimal shifts and changes in wave speed of solitary waves. We 
deduce asymptotic stability for solutions in a space of perturbations naturally constrained to omit 
these neutral-mode components, being symplectically orthogonal to them. Asymptotic stability is 
obtained in a norm that is weighted spatially to decay exponentially behind the wave profile. The 
time decay of such a norm corresponds to unidirectional scattering behavior for wave perturbations. 
The weighted-norm linear stability analysis is also used to obtain a spectral stability result in an 
unweighted energy norm. Our main results are stated precisely in section 3. 

The use of exponential weights to obtain nonlinear asymptotic orbital stability for solitary waves 
was developed for KdV solitons by Pego and Weinstein [36], for regularized long- wave equations by 
Miller and Weinstein ^31j and for Fermi-Pasta-Ulam lattice equations by Friesecke and Pego |14t [TE\ 
[TU\\n\. Finiteness of an exponentially weighted norm imposes a condition of rapid decay in front of 
the wave profile. But Mizumachi \33\ [32] recently showed how to prove asymptotic orbital stability 
for FPU solitary waves perturbed in the energy space, by using exponential weights together with 
dispersive wave propagation estimates as developed by Martel and Merle. 

Nonlinear stability for solitary water waves remains an open problem. This issue would likely 
involve a general global existence theory for small-amplitude 2D fluid motions, which is not yet 
available despite the substantial progress on well-posedness questions by Wu [121 HI] . 

There are a number of other works on (in)stability for 2D solitary water waves that concern the 
case of waves of depression with large surface tension. These include results on 2D spectral stability 
for finite- wavelength perturbations [23], spectral instability for transverse (3D) perturbations [35], 
and a full nonlinear instability result for 3D perturbations by Rousset and Tzvetkov [39]. 

A convenient tool for singular perturbation theory, used in ||35| and in the present paper to 
study spectrum in the KdV scaling limit of long time and length scales, is an operator-theoretic 
generalization of Rouche's theorem due to Gohberg and Sigal [19]. This use of the KdV scaling 
contrasts with works by Craig [12] and Schneider and Wayne [H] that concern the validity of 
the KdV approximation for water waves over time scales of order 0(e~^) for waves of amplitude 
O(e^) that are long with length scales of order e~^. Our use of the KdV approximation occurs 
in the spectral domain, where it is used to obtain partial information regarding the behavior of 
solutions to the linearized equations in the limit t — )• oo. To establish stability for time and space 
scales unrelated to the regime of validity of the KdV approximation requires a different technique 
for dealing with the linearized Euler equations, which resemble a wave equation with variable 
coefficients. We develop a method that obtains resolvent bounds from symmetrized weighted- norm 
energy estimates that use Fourier filters to cut off low frequencies. 
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2 Equations of motion and eigenvalue problem 

In this section, we derive the equations of motion hnearized around a solution steady in a frame 
moving at a constant speed c to the right, and formulate the associated eigenvalue problem. 

Basic equations. We deal with an inviscid, incompressible and irrotational fluid of constant 
density p that is bounded above by a free surface y = rj{x,t) and below by a horizontal rigid bottom 
y = —h. The velocity field (ti, v) is related to the velocity potential (f) and the stream function if) 

by 

{U,v)={(l):^,(t)y) = {ll)y,-ll)x)- (2.1) 

On the free surface y = r]{x,t), the kinematic and Bernoulli equations are: 

dtri + uria; = V, (2.2) 
dt4>+k{4'l + 'Pl)+9V = 0. (2.3) 

To make the problem non-dimensional, we let 

(x, y, t) = (hx, hy, ht/c), (77, u, v, (p, ip) = {hfj, cu, cv, chcj), chip). (2-4) 



After dropping the tildes, the equations take again the same form in the non-dimensional variables, 
with g replaced by 

7 = ^ = ^, (2-5) 

where Fr = c/ \fgh is the Froude number. 

In the fluid region, where now — 1 < y < r]{x,t), —00 < x < 00, the velocity potential and 
stream function are harmonic and are taken to satisfy the no-penetration boundary conditions 

^y{x,-l) = 0, ip{x,-l)=0 (-oo<a;<oo). (2.6) 

The dynamics is described in terms of the surface traces defined by 

^x,t) = <j){x,ri{x,t),t), ^{x,t)=il;{x,r]ix,t),t). (2.7) 



Then 



and we will write 



VJ- \^^J ytpx + Vxi'yJ \rix -^)\v 



(2.8) 



M(..)^(;-), M(..)-^^. (2.9) 

The non-dimensional equations of motion now take the form 

dtT) = v-r]^u = -V, (2.10) 
dt^ = dtcj) + (pydtrj = --i-q - \{u^ + v^) - v{v - r]xu) 

= -^ri-^{U,V)M{r^,)-\U,Vf. (2.11) 



R. L. Pego and S.-M. Sun 



5 



After transforming to a moving frame with x = x — t (= {x — ct)/h dimensionally) and dropping 
the hats, the time derivative dt is replaced hy dt — dx- A sohtary wave is a steady solution of the 
resulting equations. 

It is convenient to regard the wave motion as determined by the evolution of the pair 
with ^ and V = determined from (ry, <I>) by solving for the stream function using Laplace's 
equation and the relevant boundary conditions, namely (suppressing the t variable) 

i'xx + i^yy = {-OO < X < OO, -1 < y < r]{x)), (2.12) 

V'(x,-1) = 0, ipy - rjxi'x = U{x) {-oo < X < oo, y = r]{x)). (2.13) 

We write 

q> = Ur,<^ = ip{x, r]{x)), V = (2.14) 

Up to a normalization, "H^ is a Hilbert transform for the fluid domain. (Note + iip is an analytic 
function of x + iy.) This map will be studied in detail in a later section. 

Linearization. We linearize the equations in the moving frame about a steady solution, de- 
noting linearized variables with a dot. These linearized equations of motion take the form 

= {dt-dx)'n + v, (2.15) 

= {dt- dx)^ + -ifi + uU + vV - uvdxV- (2-16) 

Of course U = dx^- To relate V to (?),$), we linearize the boundary- value problem ()2.12p - (j2.13p 
by formally differentiating with respect to a variational parameter. The variation ip is harmonic in 
the fluid domain, zero on the bottom, and on the free surface y = r]{x) satisfies 

if = 'lpy- rjxipx - Vx'tpx + {tpyy " rix'4'xy)r]. 

Since tpyy — rjxtpxy = —dx{tpx{x,rj{x))) and —^px = v , this means 

tpy - Vxi^x = dx{^{x) - r]{x)v{x,r]{x))). (2.17) 

and by ()2.14p this means ip{x,rj{x)) = 'H'q{^ — vrj). Hence V = dx^ where 

\I' = ^p{x, T]{x))' = tp{x, r]{x)) + ipy{x, ri{x))r] = 'H^(<I> — vrf) + urj. (2.18) 

We have found it to be important (much more than merely convenient) to study the linearized 
equations of motion in terms of the combination of rj and <I> expressed as 

(j) = ^-vr]. (2.19) 

This is the surface trace of the variation of velocity potential, rather than the variation of the surface 
trace. A similar observation was made by Lannes in his treatment of well-posedness for 3D 
water waves locally in time. In terms of the pair (?), (/)), the linearized equations of motion take the 
form 
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where v' is the multipher v'{x) = dx{v{x,ri{x))). Our analysis wih show that the initial-value 
problem for the linear system (|2.20p is well-posed and (conditionally) asymptotically stable in 
a certain weighted function space. The components ?) and (j) will belong to spaces of different 
order, however, and this complicates the problem of studying stability questions directly using the 
variables (?), 4*)- 

Eigenvalue problem. Looking for solutions of ()2.20p in the form (?),</>) = e^^ {rii{x) , 4>i{x)) 
leads to the associated eigenvalue problem 

, ff\Jf!0=O- (2-21) 

7 - (1 - u)v' A - (1 - u)dx) \4>i) 

The hardest part of our analysis of the linearized dynamics involves showing that, in an appropriate 
function space, this equation has no nontrivial solutions for all nonzero A in a half plane Re A > — /3 
for some /3 > depending on the wave amplitude. 



3 Main results 

Our main result is an asymptotic linear stability result for the classical family of small-amplitude 
solitary water waves that exist for Froude number slightly more than 1, meaning 7 < 1. Asymptotic 
stability is conditional on the absence of neutral- mode components arising from translational shifts 
of the solitary wave, and wave-speed variation, as is standard. The precise results involve spaces 
with exponential weights e""^ that decay to the left (having a > 0). For a € M, we define to be 
the Hilbert space 

Ll = {f\ e-/ G l2(M)}, 

with inner product and norm 

/oo 
f{x)g{x)e^--dx, ||/||a=||e'^^/||i2. 
-00 

Also, = {/ I e"^'/ G i?''(M)} will denote a weighted Sobolev space with norm that is expressed 
in terms of the Fourier transform J-f{k) = f{k) = e~^^^ f{x) dx as 

WfWns = ||e-/||^,. = (i- I^Jl + k'r\f{k + za)\'dk 

Group velocity and weighted norms. The use of these weighted norms is motivated as 
follows. For linearization at the trivial solution 77 = $ = 0, the Hilbert transform for the fluid 
domain is the Fourier multiplier T-Iq = itanhD with symbol itanh/c (see section U]). Then the 
dispersion relation for solutions of ()2.20p with space-time dependence e*'^^-*^* is 




zu = —k lb ^7A;tanh A;. 



(3.1) 
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The solitary waves that we study travel faster than the speed of long gravity waves, meaning 
c > ^/gh and so 7 < 1. Thus, in this regime the group velocity of linear waves (relative to the 
solitary wave) is always negative: 

dw 

ik^"- 

Heuristically, linear waves scatter to the left. Our analysis makes essential use of this directionality 
by measuring perturbation size using weights e""^ with a > 0. As a simple example, the solution 
r]{x,t) = f{x + t) of the transport equation dti] = dxTj satisfies ||ry(-, i)||(j = e~*||/||a. 

In analytic terms, the isomorphism / 1— )• e"^/ from to Lp' maps a Fourier multiplier A{D) 
acting on to the weight-transformed operator e'^^ A{D)e~'^^ = A{D + ia) acting on L^. The 
L^-symbol A{k) of the former is shifted to the symbol A{k + ia) of the latter. The L^-spectrum of 
A{D) is the closure of the image of the latter symbol. This is so since the resolvent (A — A{D))~^ 
is bounded in exactly when the map /a 1— (A — A{k + ia))~^ fa{k) is bounded in , where 
fa = e""^ f- For the Fourier multipliers 



A±{D) = iD± ^/-JDt&nhD, 

which correspond to the branches of the dispersion relation ()3.ip for our water-wave problem, 
the L^-spectrum shifts from the imaginary axis into the left half-plane for small a > exactly 
because the relative group velocity is negative. The same idea underlies the use of weights to 
obtain nonlinear asymptotic stability for solitary waves of the KdV equation in [36j and of FPU 
lattice equations in [151 US Ej • 

Energy and weighted norms. Zakharov [H] showed that the water wave equations have a 
canonical Hamiltonian structure in terms of (//, $) with (nondimensional) Hamiltonian 



1 poo pTj^x) 1 poo 1 poo 

- / \V<P\^dydx + - ^rfdx = - {^{-dxnr,)^ + -irf)dx. (3.3) 
^ J —00 J —I ^ J —00 ^ J —00 



The space that we use to study asymptotic stability of the linearized system (I2.20|) is equivalent to 
a weighted linearization of this Hamiltonian about a flat surface. Namely, stability will be studied 
with (?), 4>) in the space Za = Ll>i hI'^ with norm equivalent to the norm of (?7, VI) tanh £»(/)) m 

Scaling. We study waves in the regime where the parameter 



e = (3.4) 

is small and positive. For all e in this well-studied regime, there is an even solitary-wave surface 
elevation ry with rj and surface velocity {u, v) approximately given by 

•q{x) ~ u{x) ~ €^Q{ex), v{x) ~ — e^©'(ex) 

where 

Q{x) = sech2(\/3x/2). (3.5) 
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For precise statements with estimates we use, see Theorems 15.11 and lA.ll The significance of these 
results is that we use stabihty information for the KdV sohton with the profile (j3.5p to study the 
eigenvalue problem (j2.2ip for |A| small, using the KdV scaling x = ex, A = e^A. Because of this 
scaling, we take the weighted-norm exponent to have the form a = ea, where a is required to satisfy 
< a < -v/3 to have G G L^. For convenience in analysis, our stability results are formulated with 
the tighter restriction < a < ^. The parameter a is taken as any fixed number in this range. 

Neutral modes. The solitary waves we study belong to a two-parameter family, smoothly 
parameterized by translation and Froude number (equivalently translation and wave speed c). By 
consequence, as usual the value A = is an eigenvalue of Arj with algebraic multiplicity two, with 
generalized eigenfunctions produced by differentiation with respect to x and c. Denoting these 
functions with the notation 



we have ArjZx = 0, —ArjZc = Zx- The details are developed in Appendix B. (The notation (p^ 
indicates that different choices of an integration constant are made to ensure = d~^U G .) 

Solutions of (I2.20|) that lie in the neutral-mode space spanned by Zx and Zc do not decay in 
time, naturally. A necessary condition that a solutions of (|2.20p decay in time is that it should have 
no component in this neutral-mode space. The precise spectral meaning of this (being annihilated 
by the spectral projection for the eigenvalue A = 0) can be expressed in a simple form, due to the 
canonical Hamiltonian structure of the problem. Namely, it turns out to be necessary that the 
solution be symplectically orthogonal to the neutral mode space, meaning that 



Results. Our main results concern asymptotic stability for the linearized equations in a 
weighted norm, and spectral stability in an unweighted norm. 

Theorem 3.1 (Asymptotic stability with weights) Fix a € (0, ^] and set a = ae. If e > is 
sufficiently small and rj, u, v correspond to the solitary wave profile given by Theorem \5.1\ then the 
following hold. 

(i) With domain H\ x Ha^'^ , A-q is the generator of a -semigroup in Za = L"^ x H^'^ . 
(a) Whenever ReA > —^ae^ and X 0, X is in the resolvent set of Arj- 
(Hi) The value X = is a discrete eigenvalue of A-q with algebraic multiplicity 2. 
(iv) There exist constants K > and /3 > hae^ depending on e and a, such that for all t > 0, 





(3.6) 



exp{tAr,)z\\z^ < Ke 



z\\z. 



for every initial state z = {r],(t>) that satisfies the symplectic orthogonality conditions < l^.6]) . 
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Theorem 3.2 (Spectral stability without weights) For e > sufficiently small, in the space of pairs 
{rji, 01 ) such that 



oo 

(/)i(i:>tanhi:')(;/)i +?7i < oo (3.7) 

oo 



the spectrum of the operator Arj is precisely the imaginary axis. 

The asymptotic stabiUty statement in part (iv) of Theorem 13.11 will be proved as a consequence 
of the Gearhart-Priiss spectral mapping theorem [37] by establishing that the operator Ar^ has 
uniformly bounded resolvent (A — Arj)~^ for |A| large with ReA > — |ae^, and using parts (ii) 
and (iii) to infer the resolvent restricted to the spectral complement of the generalized kernel is 
uniformly bounded for all A with ReA > — |ae^. The spectral stability result in Theorem 13.21 is 
proved using Theorem 13.11 and symmetry of the problem without weights under space and time 
reversal. 

4 Riemann mapping and Hilbert transform for the fluid domain 
4.1 Riemann stretch and strain 

We will make much use of a Riemann mapping from the fluid domain il.^ to the flat strip Qq, with 

Qri = {{x,y) : — oo < x < oo, —l<y< rj{x)}, 
^0 = {{x,y) ■ —oo < x< oo, — 1 < y < 0}. 

To denote the corresponding Riemann mapping and its inverse, we write 

(.x,y) = {Zi{x,y),Z2{x,y)), {x,y) = {zi{x,y), Z2{x,y)). (4.1) 

A key quantity is the 'Riemann stretch' at the surface, given (with its inverse h = C~^) by 

C(x) := z,{x,0), h{x) = Zi(x,r/(x)). (4.2) 

The function Z2 is harmonic in the strip VLq, with boundary conditions 

2^2(^,-1) = -!, 22(^, 0) = ??(x) := ryo^(x). (4.3) 

Taking the Fourier transform in x leads to the formula 



oo 



'^tey)-y + ^lJ'' ,2, Jj-"Ms)^s,k. (4.4) 

Using the Cauchy- Riemann equation dxZi = dyZ2, we find that the 'Riemann strain' defined by 
u = — 1 satisfies 

uj{x) = C'ix) -1 = D coth D rj{x). (4.5) 
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Integrating in x with an arbitrary constant of integration, we find that we can write 

/•oo 

({x) = X — icoth D ri{x) + Co = X— / ri{s) ds + Q 1(0)1] (x) + cq, (4-6) 

v: i 1 Qi[D) is the Fourier multipUer with symbol bounded on M given by 

^ A; cosh A; — sinh A; .,,1 , 

Qi{k) = ; — = lik"^ - coth k). 

ik smh k 

If r] is given, then since rj = rjoC,, Eq. (|4.6p is a fixed-point equation that determines C, and therefore 
h = It will turn out to be more convenient in our analysis, however, to directly study the 

Riemann strain oj, and recover other quantities such as C, and r/ from this. 

4.2 Hilbert transform 

The operator %r) admits a convenient expression in terms of the Riemann stretch C,. To see this, 
first we introduce pullback operators C# and via 

C#C/(x) = U o C(x), C*U{x) = C'C#U{x) = {U o C)(x) C'(x). (4.7) 

For later use, note that since h = C~^i the chain rule yields the simple relations 

dC# = C*d, C* = C'C# = C#(l//i'), /i*C* = C*/i* = id. (4.8) 

Write ip{x,y) = ip{x,y). Then tp is harmonic in VLq, and the boundary conditions (j2.13p transform 
to ~ ~ 

. -1) = 0, dy_±{x, 0) = Q,U{x) = dj,<^ o C)(x). (4.9) 
By Fourier transform we find 



1 ., sinhA;(w + 1) /■°° ., 

SO since ip{x^r}{x)) = 7/;(/i(x),0), after an integration by parts we find 

^'(x) = V(x,??(x)) = — / e'^^^'=\iian\ik) e-'^'^ o Q{s) ds dk. (4.10) 

J — oa J — 00 

In other words, since ^ = Hri^ we have (with denoting Fourier transform) 

-Hr, = Q^noCi^, 'Ho = itanhD = J=-^{itanhk)T. (4.11) 

Here T-Lq is the Hilbert transform for the top boundary of the strip. Though we will make no use 
of the fact, it is explicitly given in terms of an integral kernel by 

HoUix) = r koix - spis) ds, koix) = . (4.12) 

J —00 Z SlUn^TTX / Z) 
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4.3 Linearization 

To justify the later calculation of generalized eigenmodes (in Appendix B), we explain here how 
the formal linearization formula in (j2.18p follows from the representation formula in (|4.1ip . 

To proceed, start with a family of (smooth) Riemann strains to small in n and depending 
smoothly on a variational parameter, and compute 

A/ N 0-1 tanhD 

C[x)=x + o w + co, r/ = — — — Lo, r] = rioQ . 

Determine conjugate harmonic functions zi, Z2 in the strip such that (j4.3p holds and zi{x,0) = 
C{x). Then the function Z{x + iy) = zi{x, y) + iz2{x, y) yields the Riemann mapping of to 0^ 
as described above. 

Also take a family of (smooth) functions ^ (free surface velocity potential) and introduce (j) as 
the harmonic extension of ^> o into f^o satisfying dy(p = at y = —1, and -0 as the harmonic 
function conjugate to (p and satisfying = at y = —1. Then 



1 r ..fe.coshA;(y + l) , 



{x,y) = — e'"^ "f- ' I e-'^'^oC{s)dsdk. 

27r / cosh K 



oo 



oo 



Write 

T{x + iy) = (/){x,y)+iip{x,y), T = ToZ~^. 

Regarding Qq as a subset of the complex plane, we have that T and Z are analytic in and that 
T is analytic in il^ = Z{Qq). Define (p and ip to satisfy 

T(x + iy) = (p{x, y) + i'4){x, y), (x, y) G J^^. 

^{x) = 'ip{x,r]{x)) is the trace on the fluid surface and satisfies 

o C = itanhD{^ o (), 

due to \I' o (^(x) = ip{x, 0) and <I> o (^(x) = (/)(x, 0) and the boundary condition ImT = Oaty = — 1. 
Denoting the derivative with respect to the variational parameter by a dot, we have 

^' o C + C^x o C = itanhD{^ o( + (<^^ o () (4.13) 

Note that ZT' oZ is analytic in l^o and has zero imaginary part on the bottom y = —1. This means 
that the real and imaginary parts of the surface trace are related by the Hilbert transform for the 
strip. But Z = zi + iz2 = + if] on y = 0, whence (abusing notation to write ^px for ipxix^i^ix)) 
with X = C(i£)) etc.) 

^xQ + ^yV = (itanhL>)((/)^C + <Pytd- (4-14) 

Using the formulas 

fl = f]0 C + Qrja^O C, ^x = (Px + (Pyllx, "^x = ^x + IpyVx, 
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together with (j4.13p now yields 

and similarly for Combining this with (j4.14p yields 

* o C - 7)^y = itanhZ)(6 o C - ?7(/)j^), (4.15) 
and composing with h = C,~^ yields the desired linearization formula (I2.18P : 

^! - ui] = n.^{^ - vif) (4.16) 
(with {u,v) = {(pxAy) = {ipy,-ipx) ony = r]{x)). 

5 Solitary wave profiles 

In this section and Appendix A, we will give a simple self-contained account of the existence of 
small solitary waves by fundamentally the same approach as Friedrichs and Hyers |13j . establishing 
the estimates that we need regarding convergence of the scaled wave profiles in the KdV limit. 
First, note that from (j2.10p - (j2.1ip . the steady equations for a solitary wave are 

1 772 _ T/2 _|_ 2771/2 

d,v = v = d^n,^, u-^ri = -{u, v)MivrHu, vf = — ^^YT^i^ — ' ^^'^^ 

whence 

r? = ^^$, U-jr^ = ^{U^-V^)+^rjV^. (5.2) 



Using (j4.1ip and changing variables by applying (^^ we must have rj = itanh o (^), hence by 

uj = (' -1 = {D coth D)r/ = d{<!> o C) = C'U o (. (5.3) 

Then we find 

uj bp- Or] itanhL'a; 
Uo( = - = u:-- , VoC = ^ = — . 5.4 

It is convenient to apply C# to ()5.2b ). and isolate uj on the left-hand side. This turns ()5.2b ) 
into a fixed-point equation for the Riemann strain uj, in the form 

/ tanhD\"Vic<j2 - 5(itanhZ)w)2(l - 27r/)\ tanhD , , 



The following result provides scaled bounds for the fixed point approximating the sech^ KdV profile 
from (|3.5p . The proof is given in appendix A. 
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Theorem 5.1 Let a e (0, \/3), m > 2, u e (0,1), and @{x) = sech2(\/3x/2). Then for e > 
sufficiently small, equation /i5.5\) has a unique even solution in of the form 



uj{x) = e^9{€x) 

with \\9 — ©11//™ < e'^- Moreover, the map uj is smooth. 



(5.6) 



The coefficients that appear in the linearized system (j2.16p can now be expressed as follows. 
Using (|2.8p . (jS.ip and (|5.4p . on the fluid surface we have the formulas 



u 



a; + cj^ + (i tanh Duj)"^ 



l + y2 ' 



-^(1 - u) 
i + y2 ^ 

-i tanh Doj 



dC#v 



(5.7) 



(5.8) 



(1 + w)2 + (itanhZ)cj)2' (1 + a;)^ + (i tanh L»w)2 ' 1 + tj" 

Note that r/, (^u and C#^' s-^'s even functions (since functions have the same parity as their 
Fourier transform). With the choice cq = r]{s) ds, ( is odd, and rj and u are even, with 
V = (u — l)r]x odd. Formally, we have the leading order approximations 



Lo{x) ~ ri{x) ~ C#u{x) ~ e^0(e^), C#v'{x) ~ — e^0"(ex) 
For making estimates involving the quantities in (|5.8p it is useful to note that 

.tanheZ) 



|itanh eD^llji^i = ed 9 < e| 



eD 



(5.^ 



(5.10) 



//I 



6 Transforming the system 

Flattening. Given the form of Hrj in (j4.1ip . it appears convenient to transform the eigenvalue 
problem in (|2.2ip to work in variables associated with the flattened domain. We make a similarity 
transform of (j2.2ip by applying the operator = C'C# from (j4.7p to the first equation and C# to 
the second, introducing the variables 



(6.1) 



Noting that C*<9/i* = dC#h'h# = d{l/C), becomes 

-DtanhD 



7 



l — Vl 

c 



A 



1 — Ml 

c 




where (with formal leading order behavior indicated) 

ui = C#ur^ u, vi = C#-f~'^ {1 - u)v' 



0, 



;i - c#u)dc#v 

7(1+0;) 



(6.2) 



(6.3) 
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Approximate diagonalization. In order to reduce the eigenvalue problem to the 'right' scalar 
equation, it is helpful to balance off-diagonal terms (up to a commutator) and diagonalize the 
leading part of the operator. Let us define p, q, Up, Uq, and for later reference also p and Up, so 
that 

1— !^=p=l + np, = q = l + Ug, ^ = p=l + up. (6.4) 

Asymptotically we expect 

Up ~ —2a;, Uq ~ Up ~ — jo;. (6-5) 

To make precise estimates, we write 

Up{x) = e^Up{ex), Uq{x) = e^ng(ex), Up{x) = e^np(ex), 

and apply the scaled bounds from Theorem 15.11 and (jS.lOp to the expressions in (|5.8p , using 
standard calculus inequalities. Straightforward computations yield the following. 

Lemma 6.1 For e > sufficiently small, the norms ofup, Uq and Up are bounded by a constant 
K independent of e, and the functions Up, Uq, Up satisfy the pointwise bounds 

\up\ + \uq\ + \up\ < Ke"^, \up\ + \u'q\ + \u'p\ <K€^. (6.6) 

Furthermore, as e ^ we have 

||np + 2G||j^i ^ 0, WHq + ^eWni ^0. (6.7) 

Introduce the operator (Fourier multiplier) 



S = v^-7DtanhZ). (6.^ 
In order to balance orders of differentiation in the system, we change variables via 



93 



The system (j6.2p then takes the (partially symmetrized) form 

'X- dp + Ri qS \ f r]3 



(6.9) 



sq x-dp+R,)\<t.,r''' ^'-'"^ 

where Ri and R2 (which will both turn out to be negligible) are given by 

Ri = dp — qdpq~^ = (dq — qd)pq~^ = q'pq~^ , (6-11) 
R^ = dp- SpdS~^ =p' + [p,S]S-^d. (6.12) 

Finally, we approximately diagonalize by changing variables via 



1 1 / \^3j ' 2 V-1 1/ V03 



3.13) 
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Then the system (|6.10p takes the form 

X-dp \ l(-Sq-qS + Ri+R2 -Sq + qS + Ri- R2\\ (m\ ^r. 

\-dpJ ^ 2\Sq-qS + Ri- R2 Sq + qS + Ri + R2 J J ^ ' 

We make a few observations regarding the form of this system: First, the operators R2 and 
[S, q] = Sq — qS involve commutators and will turn out to be bounded, while Ri is just a multiplier. 
So the off-diagonal terms are bounded operators, involving no derivatives. Second, as is needed for 
energy estimates, we will invoke the symmetrization identity 

dp = y/p d y/p \p' . (6.15) 

The operator \{Sq + qS) can be explicitly symmetrized (for energy estimates) up to a (double) 
commutator in terms of p = yjq: 

\{Sq + qS) = y/^Sy/^+\[[S,y/^,y/^]. (6.16) 

Finally, note that the weight-transformed operator e°'^Se~°'^ is a Fourier multiplier with symbol 

S{k + ia) = V-7Ctanh^, ^ = k + ia. (6.17) 

The principal square root is used here and the real part is nonnegative. We define 

A+ = d + S, A-=d-S. (6.18) 

It is easy to see that these formulae define closed operators in with domain i/^J and with spectrum 
given by the range of the weight-transformed Fourier multipliers 



k ^ A±{i) =ii± V-7^tanh^, ^ = k + ia, /c G M. 

Final form as system. Based on these observations, it will be convenient to write the eigenvalue 
problem as follows. We use (16. 4p to write the operator in the (1,1) and (2,2) slots of (I6.14p as 
A — ^11 and A — A22 respectively, with 

Aii=A++U + Jii, U:=dup + Sug, (6.19) 

A22 = fi- + J22, B± := y/pd^ ± y/qSy/^. (6.20) 
The system (j6.14p then takes the form 

(^-^) (::)-"■ ^Kt t)' 

with the 'junk terms' Jij given in terms of Ri = q'pq^^ and R2 = p' + [p, S]S~^d by 



Ju Ji2\ ^_l(R2 + Ri + [S, q] R1-R2- [S, q] 

J21 J22 2\Ri-R2 + [S,q] Ri + \p,S]S-^d+[[S,p],p] 



3.22) 



Another way we will sometimes use to write the (1,1) component of A is 

Aii = B+ + Jii, Jii = -\{Ri + [p,S]S-^d-[[S,plp]). (6.23) 
Our main results to be proved in this paper now amount to the following. 
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Theorem 6.2 (Asymptotic stability with weights) Fix a G (0,^] and set a = ae. For e > 
sufficiently small the following hold: 

(i) With domain (H^)'^ , A is the generator of a semigroup on {L'^)'^ . 

(ii) Whenever Re A > — ^ae^ with X 0, X is in the resolvent set of A. 

(Hi) The value X = is a discrete eigenvalue of A with algebraic multiplicity 2. 

(iv) Restricted to the A-invariant spectral complement Ya of the generalized kernel of A, the semi- 
group e"^* is asymptotically stable, satisfying 

\\e-^'z\\a<Ke-^%z\\a 

for all t > and z € Ya, with some constants K > and (3 > gCte^ depending on e and a. 

Theorem 6.3 (Spectral stability without weights) For e > sufficiently small, with domain {H^Y' 
in the space (i^)^, the spectrum of the operator A is precisely the imaginary axis. 

Equivalences. The statements in these correspond directly to those in Theorem 13.11 due to the 
following facts. First, the map (?72)</'2) ^ (??4,(/'4) is clearly an isomorphism from = x h\^'^ 
to (La)^ when a > 0. (Note the symbol S[k + ia) does not vanish at A: = in this case). Second, 
the map 7]i i— t- 7]2 from (|6.ip is clearly an isomorphism on (and on H)^). For example, 

hit = r e^^lr?! oC(.)C'(.)|^e2'^(«^)-^) = {l + 0{e^))\W\\l 

J-oo C (s) 

since pointwise \C' — 1| +ae\C{s) — s\ = O(e^) uniformly, due to Theorem 15. 11 Next, the composition 
map C# is an isomorphism on H'^ for s = 0, 1 and 2, hence also for s = ^ and | by interpolation 
(see [1], particularly Theorems 3.1.2 and 6.4.4). Therefore the map 

is an isomorphism on H^- (Note that the multipliers e^'*^''*^^^"^^ = / + O(e^) on for s = 0, 1 
and 2, as is easy to check using Theorem 15.11 ) 

Finally, we claim that the transformation steps ()6.ip . ()6.9p and ()6.13p map the space of pairs 
(t/i,0i) satisfying condition (j3.7p of Theorem 13.21 isomorphicallv to the space of pairs (7/4,^4) G 

X L^. The first step to show this is to see that (174, (/)4) G x is equivalent to finiteness of 
the linearized energy: 

/oo 
M-dxnr,)(l)i + 777? dx < 00. (6.24) 
-00 

The key point here is that since (1)2 = C#<Pi and Hri = C#^^oC#5 due to ()4.8p the change of variables 
X = C{x), dx = C'{x)dx yields 

/oo foo /•oo 

M-dx'Hr^(t>idx= / {C#(t>i){-Q#d^'H^(^i)C.'dx= / </.2(Z)tanhD)</.2 = ||(/)3||i2. 
00 J — 00 J ~oo 

The second step is to demonstrate an equivalence of norms. We claim 
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Lemma 6.4 For some positive constants c_ and c+ independent of (pi, 

(/)i(L>tanhL>)(/)idx < / (t)i{-d,,'Hri)(l)i dx < c+ j (j)i{DtanhD)(j)idx. (6.25) 

-oo J —oo J —oo 

For the proof, it is enough to consider (pi smooth and rapidly decaying on M. Let cp and (p be 
functions that are harmonic in the fluid domain $7^ and the flat strip 0,q respectively, and satisfy 
the boundary conditions 

M^) = Hx,vix)) = Hx,0), = Py{x,-l) = Py{x,-l). (6.26) 

Then 

oo /-oo rT){x) 

M-dxnn)cPidx = / / \VcP\^dydx, (6.27) 

OO J — oo J — 1 

/oo /"OO /"O 

<pi{DtanhD)(pidx = / / \V^\'^dydx. (6.28) 
-oo J — oo J —1 

Moreover, the function (j) (resp. (p) minimizes the double integral in ()6.27p (resp. ()6.28p ) among 
functions satisfying the same Dirichlet boundary conditions. Let X : JIq ~^ be a smooth (but 
non-conformal) change of variables of the form X(x, y) = {x^ y{x, y)), such that X{x, 0) = {x, r/(x)). 
The function cp = p o X is smooth on and satisfies (p{x,0) = 4>{x,r]{x)) = (pi{x). Using 
the minimizing property of (p, then changing variables and using that the gradient and (inverse) 
Jacobian of X are uniformly bounded, we find 

/oo rO I'oo rO 1 roo rri{x) 

1 \V^^dydx< 1 1 \V4>\'^dydx<— / \V(p\'^dydx. 
-oo J-1 J-oo J -I C„ J -I 

This establishes the first inequality in the Lemma. The other one is similar. 



7 Estimates on commutators and junk 

In order to bound the junk terms, we need to bound commutators of S with the multipliers p, q and 
p, or equivalently with Up, Ug, and Up, since [5,1] =0. The functions tip, Ug, Up all have the scaled 
form e^G(ex), where G depends on e but remains bounded in H^. The following result provides a 
general estimate for the commutator of a Fourier multiplier and a multiplier with this scaled form. 
Write {kY = (1 + A:2)«/2, {DY = (1 + D'^f'^. 

Proposition 7.1 Let V , Q and IZ he Fourier multipliers with symbols P , Q and R respectively, 
and let s > 0. Let g{x) = e^G{ex) where G: M — )• M is smooth and exponentially decaying, and let 
/ : R — )• M 6e smooth with compact support. Then 

mQ,g]nfU2<G.CG\\f\\L^, 

where 

k.kGK {k-kY J-oo 27r 
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Proof. Using the Fourier transform and Young's inequality, since g{k) = €G{k/e), we have 

2 



P{k){Q{k) - Q(k))eG 



k-k 



R{k)fik) 



dk 
2^ 



dk 
2^ 



oo / />oo 



< c 



< C'*Cg1uiil2. 



k-k 



G 



k-k 



Corollary 7.2 Suppose < a < 7r/4. With g = Up, Ug or Up, there exists K > Q such that for all 
small enough e > 0, we have the operator norm estimates 



\\[S,g]S-^d\\a<Ke\ 

|«^ll||a + ||<^12||a + ||<^2l||a + ||'^22||a < Kc^ . 



(7.1) 
(7.2) 



Proof. Observe that for each of the indicated choices for g, we have that G is uniformly bounded 
in as a consequence of Lemma l6.ll So, using s = 2 — | we have 



Cg< 



{k)-^/^dk 



1/2 



{k)^\G{k)\^dk 



1/2 



< K 



(7.3) 



independent of e. And, the operator {D)^^^'^S^^d is uniformly bounded on since its weight- 
transformed symbol is (0~"^'^^*C/V~7C tanh^, which is uniformly bounded. Hence it suffices to 
show that with the choices P{k) = 1, Q{k) = a/— ^tanh^, R{k) = {Cl^^"^ = k + ia), we have 



To prove this estimate the idea is to show that with = k + ia, ^ = k + ia, 



\Q{k)-Q{k)\ 



^ tanh ^ — ^ tanh ^ 



Q{k) + Q{k) 



< 



K\k-k\ ^ K\k-k\ 



max 



(lJ^|l/2)- ^^)l/2 



(7.4) 



(7.5) 



and conclude through scaling by e. 

To prove (|7.5p . first note that \Q{k) — Q{k)\ < K\k — k\, since Q'{k) is uniformly bounded, as 
is easy to show. Suppose now that A; > 1, without loss. If A; < then 

\Q{k) - QCk)\ < \Qik) - Q(0)| + |Q(0) - Q(A:)| < K{\k\ + |A:|) = K\k - k\. 

If A; > 0, then one computes explicitly that 

sinh 2k + i sin 2a 



-^tanh^ = — (A + ia) 



cosh 2fe + cos 2a 
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and finds that Q{k) lies in the fourth quadrant of the complex plane together with Q{k), and so 

\Qik) + Qik)\>\Qik)\>{i)yyK, 

since tanh^ is bounded away from zero. As the map k i— t- — ^tanh^ is uniformly Lipschitz, the 
estimate (|7.5p follows. This proves the commutator estimate (jT.ip . Using this together with (j6.6p 
and the fact that Sd~^ is bounded, the remaining estimates in (j7.2p follow directly. 

8 Symbol expansions and estimates 

Here we develop basic approximations and key estimates that concern the Fourier multiplier A-\.. 

8.1 Low- frequency expansion and KdV scaling 

The Taylor expansion of tanh at zero, 

tanhe = e-k' + 0(C'), 



and the fact that for ^ with positive imaginary part one has y— ^ = yields 



V-etanhe = -i^ + liC^ + 0{f). (8.1) 
We note that if /c G M and |A;| is sufficiently small, 



< 1 - (8.2) 
In the KdV long-wave scaling, one replaces ^ in (|8.ip by e<^, and A by e^A. We find (recall 7 = 1 — e^) 



A+{eO = ieC + V-Te^tanheC = e^{^i^7i + ^ifl + fO{e^f)), (8.3) 

where 

71 = 2e-2(l - Vl-e^) = 1 + 0{e^). 
The KdV-scaled weight-transformed symbol of —A + A+ with ^ = k + ia, A = e^A = e^(Ar. + iXi) is 



€-^{-e^~X + ieC + V-7eetanhe0 = -A + i^^Ti + ^^^^7 + e^O(e^C^) 

= (-A, - ia7i + ^{a^ - 3ak'^)j) + i(-A, + iA;7i + ^(fc^ - 3^0^)7) + fO{e^f). (8.4) 

This corresponds to the purely formal KdV approximation (writing A-^ = A-^-{D)) 

-A + A+ieD) ~ e\-~X + ^8- ^8^). 
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8.2 High-frequency estimates 

Lemma 8.1 If z € C, a > 0, then Re ^/z < a if and only if ^{\z\ + Re z) < 
Proof. Write y/z = u + iv where n > 0. The result fohows from 

\z\ = \u + ivl"^ = + v'^ , Kez = u'^ — v'^. 
Lemma 8.2 Suppose ^ = k + ia with A; G M and < a < 7r/4. Then 



tanh k 



< ReV-gtanhg < — \ — - — . (8.5) 

\/cos2^V k 



Proof. By the previous lemma, if /3 > and w = y'— ^tanh^, then Re w < (3 \i and only if 

1^1 1 tanh ^1 - Re(^tanhO < 2/32. ^g g^ 

We may write 

— sinh 2k + i sin 2a u + iv 

tanh^ = — = ; = 

+ e s cosh 2k + cos 2a Dq 

with u = sinh2A;, v = sin 2a, Dq = cosh2/c + cos 2a. Then (|8.6p is equivalent to 

{k" + a2)l/2(^2 ^ ^2)1/2 _ _ < 2/32z)o, 



or (taking /c > without loss) 



2\ 1/2 / 2\ 1/2 n 



l + TI U + - -l + 7-^2/3^T^. (8.7) 
k^ ) \ ) ku ku 

Note that 

u sin 2a ^ a 
n sinh2A; ~ k 

Using this bound on the left-hand side of (j8.7p . we find that (|8.7p is implied by the bound 

Since Dq > (cosh 2k + 1) cos 2a = 2 cosh^ k cos 2a and u = 2 cosh k sinh fc, ()8.8p is implied by 

a^ tanh A: 
cos 2a k 

This yields ()8.5p as claimed. 



(8.9) 
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Corollary 8.3 Take 7 = 1 — e^, a = ea with < a < ^. For e > sufficiently small, we have that 
for all real k, with ^ = k + ia, 



Re V-7Ctanh^ < ea(l - , 



tanh A; 



(/ tanh A; \ 1 
-1 + (1 - ie^)Y — ^ < -l^^^- 

Moreover, uniformly for \ > —le^a, the operator norm of the resolvent of satisfies 



\\{X-A+rX<^. (8.10) 



The first inequality follows since 



7 1-6^ 



< — < 1 - < (1 - W 



and the resolvent bound follows since |A — ^+(^)| > Re(A — A+{S,)) > ^e^a for all k. 
For later reference, we note that for = k + ia with \a\ < vr/S, A; G M, 

, ^, I sinh 2fc| + I sin 2a| 

tanh^ < : — T < 1- (8-11 

cosh 2A; + cos 2a 

9 Semigroup generation and scalar reduction by elimination 

To start our analysis of the linearized dynamics governed by A, we use energy estimates to es- 
tablish resolvent bounds for the symmetrized operators B± that dominate the diagonal of A. By 
consequence, we show in this section that A generates a semigroup in {L'^)'^, with a = ea for 
a G [0, |]. Also we will show that if a G (0, A — A22 in (|6.2ip is uniformly invertible on for 
all A satisfying Re A > —^ea. This allows us to eliminate (/>4 in the eigenvalue problem (j6.2ip and 
reduce to a scalar, nonlinear eigenvalue equation for 774 in the form 

(A - ^11 - Ji2(A - ^22)" V2i)r?4 = 0. (9.1) 

Lemma 9.1 For some constant K independent of e, a and A, z/ a G [0, e > is sufficiently 
small, and ReA > — ea(l — Ke"^), then A is in the resolvent set of B^, with 

||(A - S_)-i|U < ^ . \ (9-2) 
Re A + ae(l — Ke^) 

and ifKeX > Kae^ then A is in the resolvent set of Bj^, with 
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Proof. The main step in the proof is to perform energy estimates for each term. Let z be 
smooth with compact support, and let Za = e°'^z. Write da = e"'^de~"'^ = d — a, and recall a = at. 
We compute (using the fact that p > 1 — Kt^ in the last step) 



oo 



R'e(-Vp5vpz,^)a = -Re/ (v^5v^z)ze^"^(ix = -Re / (^a^/v Za)^Zadx 

J — oo J — oo 

/oo 

plzal"^ dx > ae{l - Ke'^)\\z\\l. (9.4) 
-oo 

Due to (|6.17|) above, with ^ = A; + m, 

{^S^za)zae^''^ dx = Re ^/ --/Ctanh^lT^ Zaf — . (9.5) 

oo J-oo '^'^ 

By Corollary 18.31 we find < R<i{y/qS ^ z , z) a < ea\\qz\\l < ea{l + Ke'^)\\z\\l. Hence it follows that 
for all smooth z with compact support 

Re{iX-B)z,zU ^ neX + aeil-Ke% (9.6) 

\\4i 

^^((^-^^)^^^)<^ > ReA-W. (9.7) 



When the right-hand side is positive, this proves A — A± is uniformly invertible on its range, 
satisfying the respective estimates in (j9.2p and ()9.3p . 

To prove that A is in the resolvent set of B- , what remains to prove is that the range of A — B± is 
all of Lq. To accomplish this, we use a perturbation estimate to establish that a fixed value A = 1 is 
in the resolvent set for small enough e, then invoke an analytic continuation property of resolvents. 
For A = 1 > fixed, if e is small then we will show 1 — B± is a small relative perturbation of the 
Fourier multiplier 1 — A± from ()6.18p . with 

11(1 - A±r\B^ - A±)\\a < Ke^ < 1. (9.8) 

By perturbation it follows 1 is in the resolvent set of operator B± and the range of A — B± is all of 
for A = 1. Using the Neumann series for the resolvent, we see that the resolvent of any closed 
operator can be analytically continued to any set where the resolvent has a uniform apriori bound 
(see Theorem III. 6. 7 of [24J). By consequence, the resolvent set of B- (resp. B-^-) contains the 
entire right half-plane where the right-hand side of (j9.6p (resp. (|9.7p ) is positive. 
We proceed to prove (j9.8p . We compute that 

B± — A± = dup — ± Suq ± [p, S]p. 

Since Sd~^ is bounded, Corollarv 17.21 and (|6.6p imply 

\Up\ + Wq\ + \p\ + \\[S,p]p\\a < Ke^. (9.9) 
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We claim that for some constant K independent of e and a, 

11(1 - A±)-^d^\\a <K for i = and 1. (9.10) 
The symbol of the weight-transformed operator — e~"^(l — A±)~^d^ e""^ is 

m,(0 = ^ i = k + ia, (9.11) 

l-i^^y-7^tanhe 

for j = or 1. Since —i^ = —ik + a, the real part of the denominator is always greater than 1, by 
Corollary 18.31 Hence |m.o(^)| < 1 for all A; € M, so (|9.1U|) holds for j = 0. For j = 1, we have that 
|"T'i(OI — 9 ICI — 9i while for |,^| > 9 the denominator is bounded below by 

because |tanh^| < 1 by (IKTTD . Thus |mi(^)| < | for |^| > 9. Hence (f9Tn]) holds also for j = 1. 
The bound in (|9.8|) follows by combining (j9.10p with (j9.9p . This completes the proof of the Lemma. 

Proposition 9.2 For e > sufficiently small, A is the generator of a semigroup on (L^)^. 

Proof. By the Lemma just proved and the Hille-Yosida theorem, the operator 

is the generator of a semigroup on (L^)^. But ^ — is bounded, so the result follows from a 
standard perturbation theorem (see Theorem IX. 2.1 of plj). 

Lemma 9.3 For some constant K independent of e, a and A and for a G (0, if e > is 
sufficiently small then A is in the resolvent set of A22 whenever Re A > — ^ea, with 

\\{X-A22r'\\a<-. (9.12) 

ea 

Proof. Due to the bound on J22 from CorolIarv l7.2t this result follows directly from the results 
in Lemma 19.11 concerning the resolvent of B- . 



10 Resolvent bounds for |A| not too small 

For the remainder of this paper we fix a satisfying < a < ^ and write a = ea. Here we 
demonstrate a bound on the resolvent of the operator A from (j6.2ip that is uniform in A, for A in 
the right half-plane with |A| not too small. This bound, in combination with the Gearhart-Priiss 
spectral mapping theorem and our proof that the only eigenvalue of A in the right half-plane is 
A = with algebraic multiplicity 2, will allow us to obtain linear asymptotic stability in for the 
semigroup e*^*, conditional for perturbations containing no neutral-mode components. 
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Proposition 10.1 Let < < z) < 1. For e > sufficiently small, all A satisfying 

|A|>e'' and ReA > -^ae^^^'^ (10.1) 

belong to the resolvent set of A, with ||(A — < K/e^'^'^'^ , for some constant K independent 

of e and A. 

In the analysis we will make use of sharp Fourier cutoffs (Fourier filters) defined as follows. 
We specify a wavenumber threshold chosen to be k = e'', where we require v < i) < 1. Define 
projection operators tTq (low-pass), vTj (high-pass) on as follows. First, on L^, define low and 
high-pass filters by 

TTOo = J^^'^l[-f,^K]J^, TToi = I - TTqo. (10.2) 

These operators are orthogonal projections on L^. Now in L^, define orthogonal projections by 

TTo = e-""7^ooe'^^ vTi = / - vr^. (10.3) 

10.1 Resolvent bounds for for A not small 

Crucial for our estimate of (A — ^)~^ is to demonstrate uniform invertibility of the operator 

X — dp — Sq = X — All — Jii 

which appears as the dominant part of the operator in the (1, 1) slot of (j6.2ip . The aim here is to 
establish uniform invertibility of the operator above with respect to the weighted norm with weight 
a = ae. The estimate on the inverse will have the form 

\\{^-dp-SqrX<-^ (10.4) 

and be valid for A satisfying (jlO.ip . provided e > is smaller than some fixed positive constant. 
(Here and below, K is a. generic constant independent of e whose value may change from case to 
case.) Since we know ||Jii||a = 0{e^) by Corollarv 17.21 we infer that under conditions of the same 
form on A and e, A — is invertible with 

||(A-Ai)-'||a<^. (10.5) 
To prove the bound (110. 4p we study the equation 

{X-dp- Sq)z = g 

decomposing this equation in terms of Zq = tToZ, Zi = HiZ, go = nog, and gi = nig. Apply tTq 
and note noZi = 0, and the low-pass filter (nontrivially) commutes with derivatives and Fourier 
multipliers. We get 




(10.6) 
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•Aoo = A - TTo{dp + Sq)Tro, Aoi = -TTo{dUp + SUg)-Ki, 

Ao = -TTi{dUp + SUq)lTo, An = X - -Ki{dp + Sq)TTi. 

Here recall Up = p — 1, Uq = q — 1 satisfy the pointwise bounds in (|6.6|) . The low-pass Fourier filter 
satisfies \\TTod\\a < \k + ia\ < 2e^ and ||7rocS||a < Kqc^ since l^tanh^l^/^ ^ |^|| tanh^/^|^/2 < j^^^o 
for ^ = A: + ia with \k\ < k, with some constant Kq independent of e. Also ||p7ro||a + Ik^roHa < Kq. 
Now clearly, if |A| > e*^ and e is small enough so e'^ > 2Kie^ with Ki = 2Kq + Kq (we use u < u 
here), then Aoo is invertible and 

Since Sutto = (u' + 118)1:0 for u = Up and Ug, we also find (since u < I) 

WAorWa < Ke^+^, IIAolla < Ke^+^. (10.8) 

In order to establish the estimate (jlO.Sp , it suffices to show that whenever A satisfies pO.ip , An 
is invertible on with 

IIA^'IU<^ (10.9) 

This is because, after elimination of zq, (110. 6p reduces to 

(/ - A^^AioA-^Aoi)z^ = AT^\gi - AioA-}go). (10.10) 

Since 

1 1 ^ 2+' 2^ 3 

the desired estimate (jl0.4p then follows for sufficiently small e and large |A| from 

K 2K 

W^o ~'r ^ 1 1 1 1 a ~l" 1 1 1 1 a — ^1+20 y''9i\\a ~l" 1 1 1 1 a ) ^ gl+2/> 1 1 ^ 1 1 ' 

10.2 Uniform invertibility of Ai by energy estimates 

To prove the invertibility of An with the estimate (|10.9p . since z = Zq + Zi, the main step is to 
prove the energy estimate 

- Re(^i(5p + Sq)T:,z, z)a > jo^^^^^aWziWl (10.11) 

for all smooth z with compact support. Given such z, let Za = e°'^z. Recall da = e°'^de~^^ = d — a. 
Then iiQiZa is in for all m since 



[l + \k\'^T\za?dk<oo. 



\k\>k 
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1. Recall dp = ^p' and \p'\ < Ke^. Similarly as for B± which we treated before, 

/oo 
{'Ki^/pdy/pTiiz)z e^"-'^dx 
-oo 



-oo 

f CXD 



/OO 
-oo 

/oo 
plTTQiZal"^ dx > ea{l - Ke^)\\TTiz\\l. (10.12) 
-oo 



The last inequality holds since a = ea and p> I — Ke^. 
2. Now write p = y/q (as before) and compute 

TTiSqiTiZ = piTiSpzi + Cpipzi, (10.13) 

where, with Up = p — I (= 0(e^)), we can write Cpi = Cp — Cpo with 

Cp = [S, p] = SUp - UpS, Cpo = TTaSUp - UpTToS, (10.14) 

Now, since ||7ro5||a < Ke"^, evidently ||Cpo||a < Ke'^'^'^, and Corollary 17.21 implies \\Cp\\a < Ke^. 
Recall that the weight-transformed operator e"'^Se~'^^ is a Fourier multiplier with symbol given by 
()6.17p . For this we will use the high-frequency dispersion estimate in Corollary 18.31 Note that for 
> K = e'^, if e is small enough then by ()8.2p . 



>\l——>\l-^- (10-15) 

Using this with Corollary 18.31 we find 

/oo 
{piTiSpzi)zi e^"-'-^dx 
-oo 



Re / a/— 7^ isiA\^\T pzi, 

J\k\>k 



1 2 dk 

'\k\>k 27r 



> -ea{l-le^n{l + Ke^)\\z,\\l (10.16) 
3. Combining (fTOT^ with (fTOTeD and \p'\ -h \\Cpip\\a < Ke^ yields (fTOTT|) . since for small e, 
- Re{Tri{dp Sq)TTiZ, z)a 



Since An = XtTq + AuTii, it follows that if A satisfies (jlO.ip . then An has bounded inverse on its 
range with bound given by (|10.9p 

4. To prove that the range of An is all of L^, we use the same continuation approach as 
previously given for B±. We can write 

Aii = X - {d + S)7ri - TTi{dup Suq)iri. 
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For A = 1 fixed, the e-independent bound 

\\{l-{d + S)TTi)-^d7ri\\a<K (10.17) 
follows by restricting the proof of ()9.10p to frequencies |A;| > k. Then we obtain 

11(1 - {d + S)TTi)-^TTi{dUp + SUg)\\a <Ke^ <l 

for small enough e, and the invertibility of An whenever Re A > and |A| > e/K now follows as 
before for B±, by continuation based on the energy estimate in step 3. 

10.3 Bound on the resolvent of A 

Now we complete the proof of Proposition 110.11 We solve the resolvent equation 

<--^' (">(::) 

by simple elimination, writing 

04 = (A - ^22)"^(5'2 + ^21??4), 

r?4 = W,{\)-\\-Aii)-\gi + Ji2{\-A22r^g2), 

W,{\) = /-(A-^ii)-Vi2(A-^22)"V21, 

This is justified based on the estimates (|9.12p . (jlO.Sp . and the estimate 

||Ji2(A-^22)"V2i|U<i^e' (10.18) 

that follows from Corollarv 17.21 together with ()9.12p for Re A > —\iot. For the solution of the 
system, one obtains the estimates 

K / \ K / ' \ 

WViWa < -^i^ [\\9l\\a + Ke^\\g2\\aJ , \\M\a < — [\\g2\\a + K e^'"^" \\gi\\aj , 

whence the estimate ||(A — ^)^^||a < K/e^^'^'^ follows. 

11 KdV scaling and bundle limit 

It remains to study the eigenvalue problem when |A| is small, satisfying |A| < e'^. At this point we 
have shown that the eigenvalue system ()6.2ip can be reduced to the nonlinear eigenvalue equation 
()9.ip whenever Re A > —^ea. For Re A > —^e^a, we may further apply the Fourier multiplier 
(A — A+)^^ to (|9.ip . by Corollary 18.31 This reduces the eigenvalue problem to the nonlinear 
eigenvalue equation 

W{X)7]i := (/ - (A - A+r^U - J*)r/4 = 0, (11.1) 
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where 

J, = (A - A+)-\ju + Jl2(A - ^22)-V2l). 

The operator J* will be shown to be negligible. 

As we shall see, the bundle VF(A) becomes singular at A = due to the fact that zero is an 
eigenvalue of the operator A. To determine the multiplicity of this eigenvalue and establish the 
invertibility ofW{X) for nonzero A, we make use of the KdV long- wave scaling. We introduce scaled 
variables with tildes via 

x = €x, A = e^A. (11.2) 

Then d = ed, and in purely formal terms we have the following leading behavior (see (|8.3p and 
note \/—DtanhD ~ —d): 

X-A+^ e^CX - |<9 + IB^), U = dup + Sug ~ e^d{-2e) - e^d{-le). 

Thus we expect W{X) ~ VFo(A) where (with tildes omitted on derivatives) 

Wo{X) = / + (A - i5 + ^d'r^d (|G) . (11.3) 

This bundle Wo(A) is associated with the eigenvalue problem for the KdV equation scaled as 

dtf-^dj + lfdj + ldlf = o, 

linearized about the soliton profile / = = sech^(\/3x/2). 

To be clear, what we are really doing when changing variables is using a similarity transform 
in terms of the dilation operator defined by 

iTj){x) = fix/e)/V~e (11.4) 
which maps isometricaUy onto since a = ae: 

/oo roo 
|/(x)|V"^dx= / |Tj(y)|V"^dy, y = ex. 
-oo J — oo 

(Note that similarity transform does not change operator norms, but TedT~^ = ed.) 

The formal discussion above involves uncontrolled approximations in terms of derivatives. But 
this motivates the following rigorous statement in terms of convergence of bundles. Based on this 
result, the scaled operator bundle will be studied using the Gohberg-Sigal-Rouche perturbation 
theorem [lU] . 

Proposition 11.1 Define the scaled bundle W{X) = TeW{e'^X)T^^ , and let 

0, := {A e C : \e^X\ < 1, ReA > -^a}. (11.5) 
Then in operator norm on L^, we have 

sup \\WCX) -Wo{X)\\a^O ase^O. (11.6) 
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We prove this proposition by studying pieces of the scaled bundle W{X). Writing t^St^ ^ = 
^/—'yeDtan]leD = eS, t^A+t'^ = e^^+, and 

Up{x) = e^Up{ex), Uq{x) = e^Uq{ex), t^J^t^^ = J*, 

the scaled bundle is written in the form 

W{X) = / - (A - A+)-\dup + Sug) - X. (11.7) 

The proposition is implied by following convergence results in operator norm on L^, to hold as 
e — )• 0, uniformly for A € 1^^: 

WiX-A+y^B-CX-^d + ld^y^dWa^O for both ^ = a and 5, (11.8) 

||np + 2e||„ ^ 0, + ^e||„ ^ 0, (11.9) 

||J*||a^O. (11.10) 

In subsection lll.li we will prove the first limit pi.8p by studying the corresponding weight- 
transformed symbols. The third limit (jll.lU|) is treated in subsection II 1.21 The limits in (|11.9|) are 
a simple consequence of the fact that Up + 20 and itq + |0 are pointwise multipliers, so the 
operator norm is equal to the L°° norm as a function, and this is bounded by the norm, which 
tends to zero by Lemma l6.ll 

11.1 KdV limit for symbols 

Here we establish the main limit ()11.8p needed to prove the operator limit in Theorem lll.il Namely, 
we prove appropriate limits for the scaled, weight-transformed symbol of the operator A4^(X, D) = 
(A — A-^-)^^d. This symbol takes the form 

M,{X,0 = —7^--^-^^===, i = k + ia. (11.11) 
— e'^A -|- let, + V "T^C tanh e4 

The corresponding symbol for the limiting operator A^o(-^) D) = {X — i^d + ^d^)^^d is written 

-Mo(A,0 = ^-Tf-T7i- (11-12) 

-A + + giC^ 

The symbol limits that we need to prove both limits in (I11.8P are: 

|A^e(A,0 ->^o(A,OI ^ ase^O, (11.13) 



^0 ase^O. (11-14) 



These limits need to be established uniformly for ^ G M -|- and A G (The second limit will 
follow easily once we establish the first.) 
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1. First we provide simple, preliminary bounds on the limiting symbol Mq. As one sees from 
(j8.4p . for A = Xr + i\i and = k + ia, the real part of the denominator of A4q is negative, with 
magnitude bounded below by 



provided a G (0, 1) and A,. > — |a- Consequently we find 



|A^o(A,OI < 



^ eM + ia, ReA > 



-a. 



(11.15) 



(11.16) 



In particular, since \(,\ > a, the left-hand side of (|11.15p is uniformly bounded away from zero, and 
A^o is uniformly bounded. 

2. (Low frequencies) Now we carefully identify a long-wave regime where the result of Taylor 



expansion in (j8.4p yields the limits (|11.13p and (|11.14p . Fix z/q G (|, ^) and let 

Jo = e M + ia : le^l < e''"}. 



(11.17) 



Put Do = Mo{X, 0"^ E = Me{X, O'^-Dq. Then by i^, E = £,'^0{e^i'^) and since |Do| > 
we find that uniformly for ^ G /q, A G we have 



\Mo{X,i)-M,{X,i)\ 



E 



Do{Do + E) 



< 



a J l-|^|K|e^|2 



< Ke 



2uo 



(11.18) 



for small enough e, since e^l^P < e^*^" ^ = o(l). Moreover, for ^ G /q one has 



1 7 tanh 



It follows that (|11.14p holds uniformly in this regime, as well. 

3. For high frequencies the KdV limit is not relevant. In this regime, the symbols Ai^ and A^o 
must be shown separately to be small. Let us consider A^o first. When ^ G /g := M -|- \ Jq we 
have 1^1 > e''""^, and from the estimate pi.l6p it is clear that 



sup \MoCX,0\< 



a 



(11.19) 



and this tends to zero as e — )• 0. 
What remains to show is that 



sup |A(e(A,OI as e 0. 



(11.20) 



Since the square-root factor in (jll.l4p is bounded, the proof of both ()11.13p and ()11.14p will be 
complete once (|11.20p is established. This estimate is the most subtle of the symbol estimates. It 
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has nothing to do with the KdV hmit, but rather expresses a uniform stabihty property that holds 
at high frequency over moderately long time scales. Its proof breaks into two further regimes for 
1^1 and involves using Corollarv 18.31 to interpolate between low frequencies and high. 

4. (High frequencies) Here we consider the set 

Io^ = {^eW + ia:\e^\>K2}. (11.21) 
where K2 is large. In this regime, the denominator of Ai^ is estimated from below by 
|eC| - |7eCtanheC|V2 _ > l^^l - {eCl^^ - 1 > l|eC|, 

and consequently 

sup |Xe(A,e)| < 2e2. (11.22) 

5. (Transition frequencies) Fix ui with uq < vi < ^ and let 

I^={^ = k + ia: < \ek\ < K2}. (11.23) 

Recalling a = ea and 7 < 1 — o^, we apply Corollarv 18.31 together with the bound (|8.2p valid for 
\k\ small. Then we find that e small enough, with = k + ia £ Ii and —A,. < 1, the real part of 
the denominator of Aiei^^C) is negative and bounded (away from zero) by 

Re{-e^\ + zee + V-Te^tanheC) < -e^ + ea{-l + 1 - e'^"') < -^e^+^^K 

By consequence we find that 

sup |A^,(A,0| <2K2e'-'-'"S (11.24) 

and this tends to zero since vi < 

Now M + ia = /o U /i U loo, and the outstanding estimate (|11.20p is established. This finishes 
the proof of the limits (|11.13p - (|11.14p . 

Remark. The analysis of symbol limits in this section is simpler than the one carried out for 
lattice solitary waves in [17J. Partly this is due to the simple way that A appears in the denominator 
of A4f here. But partly it is due to the fact that here we must study the regime |e^A| > K hy 
other means, since the symbol A^e in (lll.lip is not bounded on the whole set where Re A > and 

11.2 Limit of junk terms 

To complete the proof of Theorem 111.11 it suffices to prove (jll.lOp . i.e., show that \\J*\\a = o{l) in 
operator norm on L^, uniformly for A G where 

= e^n, = {A G C : |A| < 1, ReA > -^ae^}. (11.25) 
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By the bound (jlO.lSp and the resolvent esthnate for ^4. in Coroharv 18.31 we have 

||(A - A+)-'Ju{X - A22r'J2i\\a < Ke\ (11.26) 

Further, Coroharies 18.31 and 17.21 implv ||(A — ^+)^^Jii||a < hence 

\\JA\a<K (11.27) 

uniformly for Re A > — ^ae^. It remains to prove that uniformly for A G Oe, 

||(A -^+)"^ Julia ^ ase^O. (11.28) 

Then ||J=i,||a = o(l) will follow. (We remark that it appears this may not be true uniformly for all 
A satisfying ReA > 0, however.) 

Because of the estimates (jll.26p and (jS.lOp and the expression for Jn in (j6.22p . it suffices to 
prove that the operators 

Jo = {\-A+)-\u'g^+[S,Uq]), (11.29) 
Ji = {\-A+r\u'p-[S,Up]S-^d), (11.30) 

are o(l) in operator norm uniformly for A G il^ — note that \q'pq~^ — y^u^l = 0{e'). 

In order to bound Jq it would suffice to note u'^ = [9, Uq] and apply Proposition 17. II with g = Uq 
and with appropriate symbols. However, the form of Ji is slightly different, so what we do instead 
is observe that the weight-transformed operators Jj = e~"'^Jje"'^ {j = 0, 1) act on a given smooth 
/ with compact support via 

{TJjf){k) = P{k)(i{k -k) + {Q{k) - Q{k))R{ky)g{k - k)f{k) — (11.31) 

for j = and 1, with 



Here ^ = k + ia, and g = Up or Uq has the form g{x) = e^G(ex) with G bounded in as in 
section [71 By almost the same short proof as that of Proposition 17.11 we find that 

II J,/||a < CM\f\\a (11.32) 

where Cq is as in Proposition 17.11 (and is uniformly bounded), and 

2\P{ek)\\ie{k - k) + {Q{ek) - Q{ek))R{eky\ 



= sup e 



sup e^\P{ek)\ 

k,keR 



i{ek — ek) 



\k-k\ 

{k-ky 



(11.33) 
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Here we take s = | as previously. This implies that the last factor in (jll.33p is bounded. To bound 
the other factors we consider the case |e^| < e'^ and its opposite, for any z> G (0, 1) fixed. 

In the first case, |e^| < e'^, the factor \P(€k)\e'^ < K uniformly since by Corollary 18.31 the real 
part of the denominator of P{ek) is bounded away from zero, satisfying ReP(eA;)^^ < — |e^a. The 
middle factor is 0{e'^) and tends to zero uniformly, since the symbols Q and R are analytic near 
^ = and Q'{k) -i and R{k) 1 as ^ = k + ia ^ 0. 

In the other case, |e^| > e'^, we note that the middle factor is uniformly bounded due to the 
estimate (|7.5|) . Due to Corollary 18.31 

ReP{ek)-' <-ReX + ea ^-l + ^^^^J < ae' - lae^+^' . 

Hence |P(eA;)|e^ < Ke'^~^'^ for small e, and we conclude that 

< i^(ei+2i> ^ 34^ 

which tends to zero uniformly for X £ 0,^. 



12 Analysis of the bundle limit 

For the remainder of the proof of our asymptotic linear stability theorem, there are two approaches 
possible. One is to proceed in a fashion similar to the treatment of FPU lattice waves in the 
KdV limit in ^7\. In that approach, one notes that any eigenfunction of A corresponding to a 
nonzero eigenvalue is orthogonal to two particular elements of the generalized kernel of the adjoint 
.4*. (In [T7| this was expressed in terms of symplectic orthogonality, using Hamiltonian structure.) 
This yields reduced orthogonality conditions that are necessary for elements of the kernel of the 
scalar bundle W{X). After an appropriate scaling, one proves convergence of these conditions to 
corresponding ones for the KdV bundle Wo{X), in a dual space. Then uniform invertibility of W{X) 
on the codimension-2 subspace satisfying the orthogonality conditions follows by a straightforward 
perturbation argument. 

We prefer to emphasize, however, that the required spectral properties follow from the bundle 
convergence theorem 111.11 bv 'soft' arguments based on the GSR perturbation theorem, and does 
not require further convergence analysis of adjoint zero modes. There are essentially only two 
'hard' points left. Namely, we need to show that (i) the bundle VF(A) is Fredholm of index zero 
for relevant values of A, and (ii) the solitary-wave degrees of freedom (translational shift and wave 
speed) naturally provide two independent elements in the generalized kernel of A. In comparing 
the need for point (i) with the alternative approach, we observe that if one knows A — ^ has empty 
kernel, one would likely need to prove a Fredholm property anyway to conclude that A — ^ is 
surjective and A is in the resolvent set of A. 

In this section we will establish point (i), and invoke Gohberg-Sigal-Rouche perturbation theory 
to characterize the null multiplicity of characteristic values of the bundle W{X). This is related to 
the algebraic multiplicity of eigenvalues of A in the following section. Point (ii) is dealt with in 
Appendix B. 
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12.1 Predholm property of the bundle 

Lemma 12.1 For e > sufficiently small, W{X) is Fredholm with index for all A G O^. 

Proof. It suffices to show that we may write 

W{X) = Wi + Wc (12.1) 

where Wi is invertible and Wc is compact. To demonstrate this, we use Fourier filters tTq and vTj 
as defined in (jl0.2p - (jl0.3p . It is convenient however to use a soft wavenumber cutoff in the range 
[k,2k] with k = where < < | (see (|11.17l) ). More precisely, fix (j){k) = 1 {\k\ < 1), 2 - \k\ 
(1 < \k\ < 2), (|A;| > 2) and set 4>,{k) = (t>{k/e'') and in place of (fTo:2]l - (fT03]l define 

7roo = -^"V.-^ TTo = e-'^^TTooe'^^ vr, = /-7ro. (12.2) 

We then define 

Wi = I -TTi,{x-A+y^u - j^, Wc = -TTo{x-A+y^u. (12.3) 

The operator Wi is uniformly invertible for A G fi^. This is so since || J*||a is uniformly small and so 
is the middle term, for we have ||9~^Z//||a < Ke^, while e~°^7ri(A — .A+)~^9e'^^' is a Fourier multiplier 
with symbol dominated by Aie, with 

\\e\.,{X-A+r^d\\a< sup |A^,(A,e)|^0 (12.4) 

as e — )• due to (|11.20p . Then, if e is small enough, Wi is invertible for all A € 0^. 

On the other hand, the operator Wc on is equivalent to the weight-transformed e~"^Wce"^ 
on L^. The latter operator is compact by the convenient compactness criterion of [34] — It is the 
sum of two terms of the form J^~^(j)iJ^(j)2, where ipi and 02 are multipliers by bounded continuous 
functions on M that approach zero at infinity. This finishes the proof of the Lemma. 

Remark. We note that in the decomposition (|12.ip . both terms Wc and Wj are analytic 
functions of A for X £ 0,^. (This fact will be used in studying the full resolvent of A.) 

12.2 Characteristic values and the Gohberg-Sigal-Rouche theorem 

We first recall some relevant basic information from [19]. (We change some terminology slightly 
for clarity. An alternative source is [18].) Let X be a Hilbert space, and suppose a function 
A I— 7- W(A) is analytic on a complex domain C C, taking values in the space of bounded linear 
operators on X, and all its values are Fredholm of index zero. A point Aq is a characteristic value 
of W if W(Ao) has a nontrivial kernel. A root vector is an analytic function z(A) with values in X 
satisfying >V(Ao)-2(Ao) = with z(Ao) ^ 0. The order of a root vector at Aq is the order of Aq as a 
zero of W(A)z(A). The null multiplicity of a characteristic value is a positive integer whose precise 
definition in general need not concern us here. The null multiplicity of Aq is always at least as large 
as the maximum order of any root vector. Furthermore, the null multiplicity equals this maximum 
order if and only if the kernel of W(Ao) is one-dimensional. 
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Suppose is a subdomain of Qq, with boundary T that is a simple closed rectifiable contour 
in r^o, and suppose W(A) is invertible for all A € F. The sum of all null multiplicities for all 
characteristic values in is denoted n(W, 0) and is called the total multiplicity of W in Q. A 
simple corollary of a far-reaching generalization of Rouche's theorem proved by Gohberg and Sigal 
[El [H] is the following. 

Theorem 12.2 Assume that for j = 1 and 2, VVj(A) is analytic and Fredholm of index zero in 
QUT. Assume that for all X £ T, VVi(A) is invertible and the operator norm 



Then W2(A) is invertible on T, and the total multiplicity n(W2,f^) = n(Wi,il). 

We apply this abstract result with Wi = Wq as defined in (|11.3p . and W2 = W as defined in 
Proposition lll.11 We take X = (i^)^, and the contour F as the boundary of the set Q = Cl^. As 
a consequence of Proposition 111.11 for e > sufficiently small, the null multiplicity n{W, 0^) = 
n(Woi^e)- But the latter number is 2, as a consequence of the following result. 

Proposition 12.3 Suppose < a < \/3 and /? = ^a{l — |a^). In L^, the only characteristic 
value ofWQ{X) with Re A > —(3 is X = 0, and this value has null multiplicity 2. 

This Proposition mainly follows from known facts concerning the eigenvalue problem for the KdV 
soliton. We provide a self-contained proof in Appendix C for the reader's convenience. 

Corollary 12.4 For e > sufficiently small, W{X) is invertible for all X on the boundary of 
17e = e^de, and the total multiplicity of W in 0^ is 2. 

13 Analysis of resolvent and eigenvalues 

It remains to complete the proof that for e > sufficiently small, the operator A has no eigenvalue 
with Re A > — |ae^ other than A = 0, which is a discrete eigenvalue with algebraic multiplicity 2. 
Conditional asymptotic stability will then follow directly from the Gearhart-Priiss theorem. 

13.1 Resolvent and spectral projection 

To begin, we note that by simple elimination, whenever A G Jl^, the resolvent equation 



||Wi(A)-i(Wi(A)-W2(A)||x < 1. 




(13.1) 



is equivalent to 



W{X)fi 

f2 



(A 
(A 



A+) \gi + JuiX - A22) ^92) 

A22y\j2lfl+92)- 



(13.2) 
(13.3) 
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Then clearly, A is in the resolvent set of A if W{X) is invertible, so any point of the spectrum of 
A inside 0^ must be a characteristic value of W{X). (The converse is not so easy to argue, since a 
root vector in need not be in H^. We finesse this point in the following argument.) 

For e > small, as a consequence of Corollary 112.41 there are at most 2 points of 0,^ in the 
spectrum of A. We claim that each such spectral point Aq is a discrete eigenvalue of A, which 
means that the associated spectral projection, 

^0 = 7^/ iX-A)-U\, (13.4) 

has finite rank. Here Tq is a small enough circle about Aq enclosing no other point of the spectrum. 
To prove this, note that from the decomposition formula ()12.ip . we can write 

W{X)-^ = W-^ - W{\)-^WcW-^, (13.5) 

from which we easily deduce from ()13.2p - ()13.3p that we can write (A — A)~^ = TZi + TZc where TZi 
is analytic in and TZc is compact. Then the integral j-^^lZidX = and it follows that Pq is 
compact. Since Pq is a projection, it has finite rank, and it follows that its range consists entirely 
of generalized eigenvectors of A. 

13.2 Algebraic multiplicity of eigenvalues 

It remains to relate the algebraic multiplicity of an eigenvalue Aq of A to the null multiplicity of 
Aq as a characteristic value of W{X). These quantities are in fact equal, but for present purposes 
it suffices to be brief and prove a simpler, weaker result. 

Proposition 13.1 For e > sufficiently small, if Aq € ile is an eigenvalue of A, then Aq is a 
characteristic value of W . Furthermore, if a Jordan chain zi, . . . ,Zk is a Jordan chain of elements 
in (H^)'^ satisfying 

{A — Xo)zj = Zj^i for j = 1, . . . , k, with zq = 0, 
then a root vector r]{X) of order at least k exists for W at Xq. 

Proof. Supposing zi, . . . , Zk is a Jordan chain for A of length k, let /(A) = Yl'j=ii^~ ^oV^^ ■ Then 
/(A) is analytic with values in {HD"^ (the domain of A) and (A — ^)/(A) = (A — Ao)'^^^ =: ^(A). 
By elimination, (|13.2p - (|13.3p hold, and consequently M^(A)/i(A) = 0(|A — Ao|^). Thus there is a 
root vector r]{X) = /i(A) of order at least k, and Aq is a characteristic value of W . 

13.3 Proof of asymptotic stability 

Recall a G (0, ^] is fixed, and take e > sufficiently small. As a consequence of the last Proposition 
and the fact from Appendix B that A = has algebraic multiplicity at least 2 for A, we conclude 
that the null multiplicity of the characteristic value A = for the bundle W is at least 2. Since 
the total multiplicity of the bundle W in 0^ is 2, we deduce that (i) there are no nonzero points of 
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the spectrum of A in Q^, and (ii) the algebraic multiphcity of A = is exactly 2. In particular, the 
kernel of A is simple (and the same is true for Vl^(O)). 

Consequently, the spectral projection Pq for A = has rank 2, and restricted to the complemen- 
tary invariant subspace Ya = {I — Pq){L'^)'^ = ker Pq, the resolvent (A — A)~^ is bounded uniformly 
for A G r^e. By consequence of Proposition llO.ll this restricted resolvent is bounded uniformly for 
all A € C with Re A > — |ae^. It follows automatically that the restricted resolvent is bounded 
uniformly in a slightly larger half-plane Re A > — /? for some /3 > |ae^. Using the Gearhart-Priiss 
asymptotic stability criterion (see Corollary 4 in [37]) gives us the conditional linear asymptotic 
stability result claimed in Theorem 16.21 

14 Spectral stability without weight 

In this section we prove Theorem 16.31 showing that in the unweighted space {L?)'^ , the spectrum 
of the operator A is the imaginary axis, if e > is sufficiently small. The proof breaks into four 
steps. For Re A > 0, we show that (i) either A is in the resolvent set or A is an eigenvalue, and (ii) 
if A is an eigenvalue in (-L^)^, then it is an eigenvalue in (L^)^- Since by Theorem 13.11 there are no 
such eigenvalues, this proves that A has no spectrum in the right half-plane. Next we show that 
(iii) A has no spectrum in the left half-plane due to a symmetry under space and time reversal. 
Finally, we show (iv) each point of the imaginary axis does belong to the spectrum of A, by a fairly 
standard construction of a sequence of approximate eigenfunctions. 

1. Suppose ReA > 0. To accomplish the first step, as in section 9 we write 

A = A,+n,, ■A* = (^Q , ^* = (j^J ■ 

By applying Lemma |9. II with a = 0, we infer that A belongs to the resolvent set of A^, and 

A - ^ = (I - 7^*(A - A^)-^){\ - A) 

We claim that 7^=f(A — ^*)^^ is compact, whence it follows that either A is in the resolvent set of 
A or it is an eigenvalue. To prove the claim, it suffices to show that each entry is a sum of terms 
each of which is a product of bounded operators, at least one of which is compact. Let L = I + d 
and note that since the domain of B± is H^^ the operators L(A — B±)~^ are bounded on L^. Thus 
it suffices to show that TZ^:L~^ is compact. By the criterion of ^3], an operator of the form gQ or 
Qg is compact on provided that 5 is a pointwise multiplier by a continuous function satisfying 
g{x) — )• as |x| — )• 00, and Q is a Fourier multiplier with continuous symbol satisfying Q{k) — )■ 
as \k\ — )• 00. 

We now deal with the various terms in TZ^^L^^ from ()6.22|) - (j6.23p . With the notation of section 
6, note q' = u'^ decays as |x| — )• 00, and has symbol (1 + ik)~^ tending to as \k\ — ?> 00. Hence 
the operator RiL~^ = pq^^q'L^^ is compact. Similarly p'L~^ is compact. 

To treat terms involving commutators, we consider first the worst term, \p,S]S~^dL~^ . Note 
that the operator S~^dL~^ is a Fourier multiplier with bounded continuous symbol ^Jk/ iSinh.k/{\+ 
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ik), hence is bounded. (The symbol decays, but we do not use this fact.) We now claim that 

[p,S] is compact. (14.1) 

We will show, in fact, that \p,S] is the uniform-norm limit of a sequence [p, 5^], where 5„ is a 
Fourier multiplier with continuous symbol of compact support. Let (/> : M — )• [0, 1] be a smooth 
cutoff function, taking the value 1 on [—1, 1] and on M \ [—2, 2] and let ip = 1 — (f>. Let Sn be the 
Fourier multiplier with symbol (f){k/n)^/—'yktanh k. Then UpSn and each compact. As in 

the proof of Corollarv l7.2l using Proposition l7.ll the operator norm of [p, S] — [p, Sn] = [up, S—Sn] 
is bounded by e^CnCc where Cq is bounded and 

fcieM {k - k)'^'^ 

Since Q is increasing and Q' decreasing for large fc, we have the uniform derivative estimate 

\{Q{k)^{k/n))'\ = \Q'{k)ij{k/n) + Q{k)^'{k/n)/n\ < \Q'{n)\ + K\Q{2n)/n\ < K/^. 

Then it easily follows C„ < Kj ^fn — )• as n — )• oo. Hence [p, 5] is compact on I?. 

Similarly the commutators and are compact, and it follows directly that TZ^L"^ is 

compact. This finishes the first step. 

2. For the second step, suppose A is an eigenvalue with ReA > and with eigenfunction 
(r/4,(/)4) G (H^)'^, the domain of A. We then can write (from (I6.14P ) 

A-^+ \ (m\ ^ (gi 

\-A-) \(t>J \g2 

with 

gi\ _ _ 1 f-dUp - SUq - UqS + Rl + R2 -SUq + UqS + Rl - R2 \ f 

92) 2\ SUq - UqS + i?l - i?2 -dUp + SUq + UqS + Rl+ R2) \(t)A 

We claim that gi and 52 lis in L\ as well as LP'. This is not difficult to check, since e~'^^Up and 
e~'^^Uq are in H^. Since ReA > 0, the Fourier multipliers (A — A±)~^ are bounded on and on 
L^. Indeed, they map the subspace L'^HL'^ of into Hiif^ (as one can check by approximation 
using smooth test functions and analyticity of the Fourier transform for < Im^ < a). It follows 
that (7/4,04) S {Ha)'^, and that A is an eigenvalue of A in the space {L"^)"^ . But there is no such 
eigenvalue for e > sufficiently small, by Theorem 13. 11 This concludes the proof of spectral stability 
for A in (L2)2. 

3. The resolvent equation for A has a symmetry under space and time reversal inherited from 
the original water wave equations. For present purposes, this is most easily studied in terms of the 
variables used in (|6.10p . for which the resolvent equation may be written (in LP x LP) 

qdpq-^ qS \ f m\ _ f fi 



(^A 2\ 

Sq x-Spds-'j v^s/ V/2; ■ ^ ' ' 
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Recall that p and q are even functions. Let C be the space reversal operator, Cf{x) = f{—x). Now, 
S preserves parity (since its symbol is even) while d reverses it. Applying space reversal to (jl4.2p 
the problem is seen to be of the same type after the replacements 

Thus A is in the resolvent set if and only if —A is. It follows at this point that the spectrum of 
A is contained in the imaginary axis. 

4. Suppose Re A = 0. Then there exists S M such that A = A+{k) = ik + i\J ^k tanh k. 

Formally, (A — A+)e'^^^ = 0. We construct a sequence of approximate eigenfunctions for A in (L^)^ 
by a cutoff and translation argument. Fix a smooth function tp with compact support, and consider 
pairs (774, (p4) of the form 

The norm of 7^4 is independent of v and r. We claim that taking = 1/n, we can choose r 
depending on n such that || (A — ^) (774, 0)11/^2 — )• as n — )■ 00. Due to the structure of A in (16.2ip it 
suffices to show that as n — )• 00. in we have (a) (A — A+)ri4^ — )• 0, and (b) {dup + Sug + Jii)??4 
and J2i??4 — ^ 0. 

To prove (a), we simply note that the Fourier transform 

^((A - A+)m)ik) = e^'^iA+ik) - A+ik))^p 

and this tends to in as — )• uniformly in r. 

To prove (b), it is convenient to note that for any fixed v > 0, — )• in as r — )• 00. 

Moreover, Upe~'^^ and Uqe~°'^ are bounded in . Then it follows, for example, that UpTji and 
UqTji — )• in as r — )• 00, and 

i?2r?4 = [diupe-"'') - 5(upe-'^^)(e''^a5L-ie-"^)) (e'^^r?4) ^ 

in L-^ as r — 7- 00, since the weight-transformed operator e"^95L~^e~"^' has bounded symbol and is 
bounded on H^. (L = 1 + 9 as above.) Similarly it follows [S, Uq]r]4 and Rirj^ — )• in as r — )• c«. 
Choosing r appropriately depending on v, this finishes the proof of (b). Thus each point of the 
imaginary axis belongs to the spectrum of A. 

A Rigorous asymptotics for solitary wave profiles 

Here we provide simple proofs of the estimates on the solitary wave profile needed for our analysis 
of the eigenvalue problem and resolvent. For a sharper treatment of solitary water waves in the 
limit e — )• see Beale's work [1]. 
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We work with the scaled form of (jS.Sp . written in terms of 9 defined through 

u{x) = e'^6{ex). (A.l) 
In terms of 9, we can write the fixed point equation (j5.5p in the form 

= F{0) = QN{9), (A.2) 
with Fourier multipHer Q and nonhnearity defined by 

, ^0'^ + e9'^ + }.(^tanheD9f(l-2je^-^^^9] 

o tanheL»\ 2 2^ M ' / , , , 

We study this equation in a weighted Sobolev space of even functions. For a > fixed, let 

= {/: M ^ R : e"^/ G H"^ , f even}, y„ = {/: R ^ M : e"^/ € i/™, / odd}, (A.4) 



with the same norm (recall (fc) = (1 + /c^)^/^) 



1/2 



One has / G (resp. Ym) if and only if / is even (resp. odd) and coshaxd^f G L^(R) for 
j = 0, ...,m. For m > 1, the space is a Banach algebra, while the bilinear product map 
if 7 9) fo is continuous from Ym x Ym to Xm- The intersection of exponentially weighted 
spaces is the direct sum of Xm and Ym'- 

n =Xm(B Ym. 

Due to the Taylor expansion of tanh, the symbol of Q has the expansion 

Q(0 = ^ , $ = k + ia. (A.5) 

Formally, the limit of the fixed point equation (1A.2P is 

9 = QoNo{9), QQ = (i-ld^y\ No{9) = l9\ (A.6) 

Provided < a < \/3, this fixed-point equation is satisfied by the KdV traveling-wave profile 

e(x) = sech2(^/3x/2). (A.7) 

This fixed point is nondegenerate in the space Xm- Indeed, the linearized map 9^9 — (5o(300) 
has bounded inverse on Xm, for the following reason. It is straightforward to show that the map 
QqQ is compact on H"^ n H^,^ (using [34]). So if ^ 1— )• ^ — Qq{?>Q9) is not an isomorphism on 
Xm-, then it vanishes for some nontrivial 9. By a simple bootstrapping argument, this 9 must be a 
smooth function satisfying (/ — + 30)0 = with e"^'0 G L^. But only constant multiples of the 
odd function 9 = Q' have this property; from standard results for asymptotic behavior in ordinary 
differential equations, any independent solution 9 grows like e^'^' as x — )• ±00. 
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Theorem A.l Fix m > 2, a G (0, \/3), CLnd v G (0, 1). Then for e > sufficiently small, equation 
iA.^) has a unique fixed point that belongs to and satisfies \\0 — Q\\xm < e*^- This fixed point 6 
depends smoothly on e. 

To prove this result, we wiU invoke a standard fixed-point lemma in the simple quantitative 
form from the appendix of p!^. To make the estimates needed, we single out one difficult nonlinear 
term and let 

N i[e) = {t&nh. eDOf. (A.8) 
Note ianheDO is odd if 9 is even, and tanheZ) maps Xm to Ym continuously. Then we write 

Q = Qo + Qi, iV = iVo + iVi + iV2, F = Fo + Fi+F2 + F3, (A.9) 

with 

Fo = QoiVo, Fi=QoA^i, = QoiV2, F^ = QiN. 

We will prove that for each j = 0,1,2,3, Fj is a smooth map on X^-, and will prove that for 
5 = e'^ > d small {u £ (0, 1) fixed) and Bs a (5-ball about G in Xm, 

\\Fjm\x^ < S, sup \\F;m\c(x^) <S, j = 1,2,3. (A.IO) 

oeBs 

(Here || • denotes the operator norm on X^.) The estimates in Theorem I A . 1 1 follow directly 

from Lemma A.l of by these estimates and the fact that Fq is smooth and / — Fq{Q) has 
bounded inverse. 

We proceed to prove the estimates in (jA.lOp . It is clear that each Nj (j = 0, 1, 2, 3) is smooth 
in Bg. Also it is not hard to see that for some constant independent of e, the remainder term in 
the nonlinearity satisfies 

||iV2(e)||x„ + sup mmWcix^) < Ke (All) 
e&Bs 

Since Qq is bounded on Xm, the estimates (jA.lOp hold for F2 with 6 = Ke. We now need one more 
symbol estimate. 

Lemma A. 2 Let u G (0, |). Then there exists K such that for e sufficiently small, 

\\Qi\\c(x^)=sup\Qi{k + ia)\<Ke^'^. 

km 

Proof. We fix u G (0, ^) and consider separately the low- frequency case |e^| < 3e'^ and its high- 
frequency complement, with ^ = /c + ia. 

1. Consider first the regime |eC| < Se'^. Let 

^0 = Qoicr' = 1 + 5^', = Qoicr' - Qicr' = eo{^e), 

so that 
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From the Taylor expansion used in (jA.SP we infer that 

E 



2!y|c|2 

\Qim<K 



Dn 



\Dr)\ 



2. In the regime |e^| > 3e'^, we estimate Qo and Q separately. First, for e < 1 we have > 9 and 
consequently |Qo(OI is handled by the estimate 

6 

It remains to bound |Q(^)|- We calculate that for = k + ia with e small, 
tanhe^ A; sinh 2e/c + a sin 2ea e^^ 



|Qo(6l<^<262-2-. (A.12) 



Re 



cosh 2ek + cos 2ea k'^ + 
sec ea sinh 2ek 2c? 
- ek cosh 2ek + 1 ^ |^ 



tanh ek 2c? 
7^ + ^ 



sec ea "^";7'" + ^ < (1 + e2a2)(i _ ^2^) + g2-2,._ ^^^^g^ 

This implies that 



ReQ(0"' > - (1 + e2a2)(i _ _ ^2-2-) > ie2--2 

for small enough e (since 2v <2 — 2v\ and hence 

\Q{i)\ < 2e2-2-. (A.14) 

This finishes the proof of the Lemma. 

Prom this Lemma, the estimates (jA.lOp for F3 clearly follow with 6 = Ke^^ . 

It remains to prove (jA.lOp for Fi, with 5 = Ke^. To do this it is convenient to note that the 
operator Qq gains regularity — it is a bounded map from Xm-i to Xm- Since A''i is quadratic, it 
then suffices to prove that for some constant K independent of e, we have 

||(tanheZ)0i)(tanhei?02)||x„_i < Ke^OiWxjmXr. (A.15) 

for all 9i, 02 £ Xm- But this follows easily since the bilinear product map is continuous from 
Yrn-i X Yjn-i to Xm-1, and 

I tanh e(A; + ia)| | tanh e(A; + 2q)| 
sup — < e sup — — < Ke (A. 16) 

km W km {(-k) 

which implies that for all 6 S Xm, 

\\{i^nheDe)\\y^_,<Ke\\e\\x^. 

This finishes the proof of the estimates in (|A.10p . 

That the fixed point is a smooth function of e is a standard consequence of the easily verified 
fact that the map (e, 9) 1— )• F{9) is smooth. 



R. L. Pego and S.-M. Sun 



43 



B Neutral modes and adjoints 

Here we verify that the translational and wave-speed sohtary-wave degrees of freedom naturally 
yield two independent elements of the generalized kernel of ^ in L^, and we demonstrate that the 
symplectic orthogonality conditions ()3.6p transform precisely to the condition that the initial data 
for the linearized equations lie in the spectral complement to this generalized kernel. 

1. Recall that by Theorem 15.11 we have a smooth family of solitary waves {rj, U) that are 
solutions of the equations (|5.ip . We exploit invariance with respect to translation by differentiating 
in X to obtain an eigenfunction of (j2.2ip corresponding to A = 0. This is slightly tricky due to the 
meaning of the variable (pi in (12.2ip . We claim that the eigenfunction has the form 

where (px = — vVx = u is evaluated on the surface (x,7/(x)). To justify this statement, we note 
that 

(1 - u)r]x = -HriU = -V. (B.2) 

To see this, recall from (|5.2p that ij = ^ = Hri^, and this equation continues to hold for translated 
wave profiles. Differentiating with respect to the translation parameter we have rj = r]x = ^ and 
^ = = (px + vr]x- Then ()B.2p follows from the linearization formula ()2.18p . 

Using ()B.2p together with direct differentiation of (jS.ip (as in (j2.20p . noting V = r]x) yields 

-a(l-n) f^^. )M=o. (B.3) 
7-(l-ii)u' -{I - u)dx J \(t)x J 

Thus Ar^zi = 0. Carrying out the transformations (16. Ij) . (16. 9p . (16.13P that lead to (I6.2ip . we let 

Due to the regularity from Theorem I A . 1 1 and the formulae (|5.8p and (|6.4p . 24 G (H^ n H\^)'^ with 

Azi = 0. 

2. Next we exploit wave-speed variation to find a generalized eigenfunction in for A = 0. 
To calculate this, it is convenient to unscale the wave speed and keep 7 at a fixed value 7 when 
computing variations. For some fixed c > ^/gh set 

qH c / 7 

7 = ^, c=- = W-, ri^{x;c) = r]{x;-f), U^{x;c) = cU{x;-f), (B.5) 
c y 7 

<^+{x;c) = c<l>+{x;j)=d-^U^=c uj{s;-f) ds. (B.6) 

J +00 



44 



Asymptotic linear stability of solitary water waves 



Then the unsealed sohtary-wave profile (t/*,^'^) is a smooth function of c that takes values in 
X Ha^'^ and satisfies 

-CTj^+Tir^^^t = 0, -cU,+^rj, + ^{U,,V,)M{7],)-\U,,V,f = 0, (B.7) 

where 14 = dx7irj,^t = cdxT]^. For c = 1 this simply corresponds to (|5.ip - ()5.2p with 7 = 7. 

We differentiate with respect to c, then set c = 1 and drop the star subscripts. Denoting the 
c-derivative by the subscript c, using the linearization formula ()2.18p we find Vc = dxJ} + dxJic and 

-'i]-r]c + Ur^i^t - vqc) + uric = 0, 
-U -Uc + JVc + uUc + v{dx'n + dxVc) - uvdxTjc = 0. 
Since U — rjxV = u = (px, this yields 



-dxil - u) OxUr, \ [ Vc\ ^ fVx 

j-{i-u)v' -{i-u)dxj ycptj ~ 



With 



c / 

rl\2 



we have — ^r^yi = zi, and find that 7/4 G (H^^) with 

-^?/4 = ^4. (B.IO) 

Adjoint modes. It is a standard fact of operator theory that the space Ya, the kernel of the 
spectral projection Pq in (|13.4p . is the subspace annihilated by the generalized kernel of the adjoint 
A* . This generalized kernel is two-dimensional (since the generalized kernel of A is) , and we aim to 
show that the annihilation conditions correspond to the symplectic orthogonality conditions ()3.6p . 

We will work with the Banach space dual of L^, and note that for the Fourier multiplier 
S = \J —^D tanh £), the adjoint is given formally hy S* = S acting in IPl^- To see this, take smooth 
test functions / and g and write fa = e"^/) 9-a = and 5a = e"'^Se~°'^. Then since the 

symbol of S satisfies S{k + ia) = S{k — ia) we have 

/oo j-OO POO ^ ^jjj. 

{Sf)gdx = I Safa¥^dx= I S{k + ia)fa{k)g^a{k) — 
-00 J —00 J —00 



00 J —00 

00 /•oo 



faS-aQ-adx = I fSgdx. 

^ — 00 

To describe the generalized kernel of the adjoint A*, first note that with the definition 

^C-Hx) 



(x;c) = c$ (x;7) = c / uj{s;j)ds, (B.ll) 

J —00 



we can repeat the arguments leading up to (jB.SP with $ and Hta replacing and H^. Then 
(jB.SP holds with = $~ — vr]c replacing 
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Next, it is convenient to work with the variables used in (j6.9p - (|6.10p and note that 

1)2' -^^'ySq SpdS-^ 

The operator ^3 acts in (L^)^- can be found by transformation from the original canonical 
Hamiltonian structure, this operator admits the factorization 

^ _ ( 1 -q-^pdS-^ 



The adjoints are given by J* = — J , C* = £, and ^3 = —CJ, acting in Then with 

one can check directly that Zg, yg G and 

* _ f-q~^pdq -qS \ f -{-fq^^Ci^^, 
^ ^ ~ V -S-^dpS) V 5-iC*r/. 

(79)-^ \ / -pd ^q^\ /C#0: 
\^L>tanhL> dp ) ^C??, 

and similarly —A^y\ = z^. Thus the generalized kernel of ^3 is the span of and 2/3. 

Now, corresponding to an arbitrary element ii = (?), 0) G = L'^xHa^'^ is is = {'jqC^r], 5C#0) S 
(L^)^. Then the conditions that and annihilate is transform as follows: 

/oo poo 
ilQC*'n)ilQy^C#4>x - iSC#<j))S-'^C*Vx dx= r]4>cc - <j)r]^ dx, (B.15) 
-00 J —00 

/oo roo 
{iqC*r]){'yq)~'^C#(l)c - {SC#<j))S-'^C*r]cdx = V<Pc - hcdx. (B.16) 
-00 J ~oo 



This shows that the symplectic orthogonality conditions (j3.6p transform to the precise condition 
that the initial data for the semigroup e"^* lie in the space Ya = ker Pq that is the spectral comple- 
ment of the generalized kernel of A. 

C Characteristic values for the KdV bundle 

Here we provide the proof of Proposition 112.31 concerning the characteristic values of the bundle 
Wo{\) = I+{\-\d + \d^y^d{lQ), G = sech2(V3x/2). 
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(In this section we will drop the tilde on A for convenience.) For Re A > (3 = \ol{\ — |a^), the 
weight-transformed operator 

e--(i _ T^o(A))e-"" = (-A + i(9 - q) + i(5 - af)-\d - q)(|G) 

is compact on and has range in H^, due to estimates similar to (jll.l5p - (|11.16p . On L^, therefore, 
Wo (A) is Predholm of index zero. If Wo{X)f = for some nonzero f ^ L'^, then e"^f E i/™ for all 
m by an easy bootstrapping argument, and / satisfies the ordinary differential equation 

(^x-id+'^d^)f + d{lef)=o. (c.i) 

By standard results, such an equation has a solution / ~ e'^^ as x — )• oo for each /i that satisfies 

A - i/x + i^u^ = 0. (C.2) 

For A > large, this equation has one root with Re/i < —a and two with /i with Re /i > —a. With 
/i = —a + it (t G M), the curve 

II- = -a+ - at^ + i{t + - a^t) 

has increasing imaginary part and real part less than —a + = —2/3. For Re A > — /3, then, 
()C.2p has a unique and simple root satisfying Re ;U < —a, hence (jC.ip has a unique solution (up 
to a constant factor) satisfying e°^/ — )• as x — )• oo. In particular, one may check explicitly (and 
easily by computer) that 

/ = a, (e''"((V3 + /.)2-(^/3 + ^ + /ie^")^/3sech2(^/3x/2))) . (C.3) 

Since Re// < —a, clearly e°^^/ G is impossible unless \/3 + ^ = 0, meaning A = and / = dxQ- 
From the analysis so far, we see that the kernel of VFo(O) in -^q is one-dimensional. To finish 
the proof, we need to show that there is a root vector at with order 2, and no root vector of order 
3. Any root vector /(A) at A = may be taken in the form /(A) = /o + A/i -|- A^/2 + O(A^) with 
/o = d^Q. And Wo(A) = Wq + VF^A + \W^\^ + 0{\^) where Wq = VFo(O) and 

wl, = -{-\d + id'y^dil®), h< = + Id'r'dilQ). 

To find a root vector of order 2, it suffices to find /i G such that Wq/i + W^o/o = 0. Since 
Wo/o = we have 

so it suffices to find /i such that 

(-ia + i93)/i + a(|G/i) + /o = o. 

Such a function can be found by differentiating the equation satisfied by the KdV wave profile with 
respect to wave speed. The function (pb{x) = 5sech^ ^/3bx/2 satisfies (pi = Q and 

i-^bd, + ldl)iPb + d,i^^l) = o. 
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Differentiating with respect to b and setting b = I, we find that 

+ + ld,e)dbipi - \d^^i = o. 

From this we see that the choice /i = —■^db(pi works and yields a root vector of order 2. This 
choice is unique up to adding a scalar multiple of /q. 

To show that there is no root vector of order greater than 2, it suffices to show that with /i as 
above, there is no /2 G such that 

Wof2 + W^i + iiy^'/o = 0. (C.4) 

If such an /2 exists, then a bootstrapping argument involving the decay estimate ()11.16p shows 
that e"^/2 G for all m. Because of the equations satisfied by /i and /o, we find that 

Therefore /2 must be a smooth solution of 

(-^5x + |93 + |a,e)/2 + /i = o. (C.5) 

Now, the function ipi = Q has e~°'^ipi in H"^ for all m and satisfies 

(-ia, + + led^)ipi = 0. 

Multiplying (jC.Sp by and integrating by parts, we find that the terms involving /2 vanish. Thus, 
for /2 to exist, it is necessary that ipid^ipi = 0. But 

/_ ^^^^^^ = 1 /_ = 1^'^' /_ > 

Hence, /2 cannot exist as required, and this proves that the characteristic value A = has null 
multiplicity 2. 
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